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Preface 

Electromagnetism is fundamental to the whole of electrical and electronic engineering. It provides the 
basis for understanding the uses of electricity and for the design of the whole spectrum of devices 
from the largest turbo-alternators to the smallest microcircuits. This subject is a vital part of the 
education of electronic engineers. Without it they are limited to understanding electronic circuits in 
terms of the idealizations of circuit theory. 

The book is, first and foremost, about electromagnetism, and any book which covers this subject must 
deal with its various laws. But you can choose different ways of entering its description and still, in 
the end, cover the same ground. I have chosen a conventional sequence of presentation, beginning 
with electrostatics, then moving to current electricity, the magnetic effects of currents, 
electromagnetic induction and electromagnetic waves. This seems to me to be the most logical 
approach. 

Authors differ in the significance they ascribe to the four field vectors E, D, B and H. I find it 
simplest to regard E and B as 'physical' quantities because they are directly related to forces on 
electric charges, and D and H as useful inventions which make it easier to solve problems involving 
material media. For this reason the introduction of D and H is deferred until the points at which they 
are needed for this purpose. 

Secondly, this is a book for those whose main interest is in electronics. The restricted space available 
meant that decisions had to be taken about what to include or omit. Where topics, such as the force on 
a charged particle moving in vacuum or an iron surface in a magnetic field, have been omitted, it is 
because they are of marginal importance for most electronic engineers. I have also omitted the chapter 
on radio-frequency interference which appeared in the second edition despite its practical importance. 

Thirdly, I have written a book for engineers. On the whole engineers take the laws of physics as 
given. Their task is to apply them to the practical problems they meet in their work. For this reason I 
have chosen to introduce the laws with demonstrations of plausibility rather than formal proofs. It 
seems to me that engineers understand things best from practical examples rather than abstract 
mathematics. I have found from experience that few textbooks on electromagnetism are much help 
when it comes to applying the subject, so here I have tried to make good that deficiency both by 
emphasizing the strategies of problem-solving and the range of techniques available. A companion 
volume is planned to provide worked examples. 

Most university engineering students already have some familiarity with the fundamentals of 
electricity and magnetism from their school physics courses. This book is designed to build on that 
foundation by providing a systematic treatment of a subject which may previously have been 
encountered as a set of experimental phenomena with no clear links between them. Those who have 
not studied the subject before, or who feel a need to revise the basic ideas, should consult the 
elementary texts listed in the Bibliography. 
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The mathematical techniques used in this book are all covered either at A-level or during the first year 
at university. They include calculus, coordinate geometry and vector algebra, including the use of dot 
and cross products. Vector notation makes it possible to state the laws of electromagnetism in concise 
general forms. This advantage seems to me to outweigh the possible disadvantage of its relative 
unfamiliarity. I have introduced the notation of vector calculus in order to provide students with a 
basis for understanding more advanced texts which deal with electromagnetic waves. No attempt is 
made here to apply the methods of vector calculus because the emphasis is on practical problem- 
solving and acquiring insight and not on the application of advanced mathematics. 

I am indebted for my understanding of this subject to many people, teachers, authors and colleagues, 
but I feel a particular debt to my father who taught me the value of thinking about problems 'from 
first principles'. His own book. The Electromagnetic Field in its Engineering Aspects (2nd edn, 
Longman, 1967) is a much more profound treatment than I have been able to attempt, and is well 
worth consulting. 

I should like to record my gratitude to my editors. Professors Bloodworth and Dorey, of the white and 
red roses, to Tony Compton and my colleague David Bradley, all of whom read the draft of the first 
edition and offered many helpful suggestions. I am also indebted to Professor Freeman of Imperial 
College and Professor Sykulski of the University of Southampton for pointing out mistakes in my 
discussion of energy methods in the first edition. 

Finally, I now realize why authors acknowledge the support and forbearance of their wives and 
families through the months of burning the midnight oil, and I am most happy to acknowledge my 
debt there also. 

Richard Carter 
Lancaster 2009 
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1. Electrostatics in free space 

Objectives 

• To show how the idea of the electric field is based on the inverse square law of force between 
two electric charges. 

• To explain the principle of superposition and the circumstances in which it can be applied. 

• To explain the concept of the flux of an electric field. 

• To introduce Gauss' theorem and to show how it can be applied to those cases where the 
symmetry of the problem makes it possible. 

• To derive the differential form of Gauss' theorem. 

• To introduce the concept of electrostatic potential difference and to show how to calculate it 
from a given electric field distribution. 

• To explain the idea of the gradient of the potential and to show how it can be used to calculate 
the electric field from a given potential distribution. 

• To show how simple problems involving electrodes with applied potentials can be solved 
using Gauss' theorem, the principle of superposition and the method of images. 

• To introduce the Laplace and Poisson equations. 

• To show how the finite difference method can be used to find the solution to Laplace's 
equation for simple two-dimensional problems. 
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1.1 The inverse square law of force between two electric charges 

The idea that electric charges exert forces on each other needs no introduction to anyone who has ever 
drawn a comb through his or her hair and used it to pick up small pieces of paper. The existence of 
electric charges and of the forces between them underlies every kind of electrical or electronic device. 
For the present we shall concentrate on the forces between charges which are at rest and on the force 
exerted on a moving charge by other charges which are at rest. The question of the forces between 
moving charges, which is a little more difficult, is dealt with in Chapter 4. 

The science of phenomena involving stationary electric charges, known as electrostatics, finds many 
applications in electronics, including the calculation of capacitance and the theory of every type of 
active electronic device. Electrostatic phenomena are put to work in electrostatic copiers and paint 
sprays. They can also be a considerable nuisance, leading to explosions in oil tankers and the need for 
special precautions when handling metal-oxide semiconductor integrated circuits. 

The starting point for the discussion of electrostatics is the experimentally determined law of force 
between two concentrated charges. This law, first established by Coulomb (1785), is that the force is 
proportional to the product of the magnitudes of the charges and inversely proportional to the square 
of the distance between them. In the shorthand of mathematics the law may be written 

Y= Q1Q2 ^ ^^^^ 

where Q\ and Q2 are the magnitudes of the two charges and r is 
the distance between them, as shown in the figure on the right. 
Now force is a vector quantity, so Equation (1.1) includes the 
unit vector r which is directed from Qx towards Q2 and the 
equation gives the force exerted on Q2 by Q\. The force exerted 
on Qx by Q2 is equal and opposite, as required by Newton's 
third law of motion. 

Examination of Equation (1.1) shows that it includes the effect of the polarity of the charges correctly, 
so that like charges repel each other while unlike charges attract. The symbol denotes the primary 
electric constant; its value depends upon the system of units being used. In this book SI units are 
used throughout, as is now the almost universal practice of engineers. In this system of units is 
measured in Farads per metre, and its experimental value is 8,854 x 10'^^ F m"^; the SI unit of charge is 
the coulomb (C). 
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Electric charge on a macroscopic scale is the result of the accumulation of large numbers of atomic 
charges each having magnitude 1.602 x 10'^^ C. These charges may be positive or negative, protons 
being positively charged and electrons negatively. In nearly all problems in electronics the electrons 
are movable charges while the protons remain fixed in the crystal lattices of solid conductors or 
insulators. The exceptions to this occur in conduction in liquids and gases, where positive ions may 
contribute to the electric current. 

1.2 The electric field 

Although Equation (1.1) is fundamental to the theory of electrostatics it is seldom, if ever, used 
directly. The reason for this is that we are usually interested in effects involving large numbers of 
charges, so that the use of Equation (1.1) would require some sort of summation over the (vector) 
forces on a charge produced by every other charge. This is not normally easy to do and, as we shall 
see later, the distribution of charges is not always known, though it can be calculated if necessary. 
Equation (1.1) can be divided into two parts by the introduction of a new vector E, so that 

E = --^r (1.2) 

and F = Q^E (1.3) 

The vector E is known as the electric field, and is measured in volts per metre in SI units. The step of 
introducing E is important because it separates the source of the electric force (Qi) from its effect on 
the charge Q2. The question of whether the electric field has a real existence or not is one which we 
can leave to the philosophers of science; its importance to engineers is that it is an effective tool for 
solving problems. 

The electric field is often represented by diagrams like Fig. 1.1 in which the lines, referred to as 'lines 
of force', show the direction of E. The arrowheads show the direction of the force which would act on 
a positive charge placed in the field. The spacing of the lines of force is close where the field is strong 
and wider where it is weak. This kind of diagram is a useful aid to thought about electric fields, so it is 
well worth while becoming proficient in sketching the field patterns associated with different 
arrangements of charges. We shall return to this point later, when discussing electric fields in the 
presence of conducting materials. 
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Fig. 1.1 The electric field of a point charge can be represented diagrammatically by lines of force. The 
figure should really be three-dinnensional, with the lines distributed evenly in all directions. 

In order to move from the idea of the force acting between two point charges to that acting on a charge due to 
a whole assembly of other charges it is necessary to invoke the principle of superposition. This principle 
applies to any linear system, that is, one in which the response of the system is directly proportional to the 
stimulus producing it. The principle states that the response of the system to a set of stimuli applied 
simultaneously is equal to the sum of the responses produced when the stimuli are applied separately. 

Equation (1.2) shows that the electric field in the absence of material media ('in free space') is proportional 
to the charge producing it, so the field produced by an assembly of charges is the vector sum of the fields 
due to the individual charges. The principle of superposition is very valuable because it allows us to tackle 
complicated problems by treating them as the sums of simpler problems. It is important to remember that 
the principle can be applied only to linear systems. The response of some materials to electric fields is non- 
linear and the use of the principle is not valid in problems involving them. 
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Before discussing ways of calculating the electric field it is worth noting why we might wish to do it. 
The information might be needed to calculate: 

• the forces on charges; 

• the conditions under which voltage breakdown might occur; 

• capacitance; 

• the electrostatic forces on material media. 

The last of these is put to use in electrostatic loudspeakers, copiers, ink-jet printers and paint sprays 
and is of growing importance in micro-mechanical devices. 

1.3 Gauss' theorem 

Figure 1 . 1 shows electric field lines radiating from a charge in much the same way that flow lines in 
an incompressible fluid radiate from a source (such as the end of a thin pipe) immersed in a large 
volume of fluid (shown in Fig. 1.2). Now in the fluid the volume flow rate across a control surface 
such as S, which encloses the end of the pipe, must be independent of the surface chosen and equal to 
the flow rate down the pipe, that is, to the strength of the source. 




Fig. 1.2 The flow of an incompressible fluid from the end of a thin pipe is analogous to the electric 
field of a point charge. It is necessary for the end of the pipe to be well away from the surface of the 
fluid and the walls of the containing vessel 

To apply this idea to the electric field it is necessary to define the 
equivalent of the flow rate which is known as the electric flux. The figure 
on the right shows a small element of surface of area dA and the local 
direction of the electric field E. The flux of E through dA is defined as the 
product of the area with the normal component of E. This can be written 
very neatly using vector notation by defining a vector dA normal to the 
surface element. The flux of E through dA is then just E-dA ^ E dA cos 9. 
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Now consider the total flux coming from a point charge. The simplest choice of control surface 
(usually called a Gaussian surface in this context) is a sphere concentric with the charge. Equation 
(1.2) shows that E is always normal to the surface of the sphere and its magnitude is constant there. 
This makes the calculation of the flux of E out of the sphere easy - it is just the product of the 
magnitude of E with the surface area of the sphere: 



Thus the flux of E out of the sphere is independent of the radius of the sphere and depends only on the 
charge enclosed within it. It can be shown that this result is true for any shape of surface and, by using 
the principle of superposition, for any grouping of charges enclosed. The result may be stated in 
words: 

The flux ofE out of any closed surface in free space is equal to the charge enclosed by the surface 
divided by sq- 

This is known as Gauss' theorem. It can also be written, using the notation of mathematics, as 



where 5* is a closed surface enclosing the volume Fand p is the charge density within it. Equation 
(1.5) looks fearsome but, in fact, it is possible to apply it directly only in three cases whose symmetry 
allows the integrals to be evaluated. Those cases are: 

• Parallel planes 

• Concentric cylinders 

• Concentric spheres 

To show how this is done let us consider the case of a long, straight, rod of radius a carrying a 
uniform charge q per unit length. From the symmetry of the problem we can assume that E is 
everywhere directed radially outwards and that the magnitude of E depends only on the distance from 
the axis. This is not valid near the ends of the rod but the problem can be solved in this way only if 
this assumption is made. The next step is to define the Gaussian surface to be used. This is chosen to 
be a cylinder of radius r and unit length concentric with the charged cylinder with ends which are flat 
and perpendicular to the axis as shown in Fig. 1.3 



flux of E = 




(1.4) 



(1.5) 



S V 
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Fig. 1.3 A Gaussian surface for calculating the electric field strength around a charged rod. 

On the curved part of the Gaussian surface E has constant magnitude and is everywhere perpendicular 
to the surface. The flux of E out of this part of the surface is therefore equal to the product of E and 
the area of the curved surface. On the ends of the cylinder E is not constant but, since it is always 
parallel to the surface, the flux of E out of the ends of the cylinder is zero. Finally, since the Gaussian 
surface is of unit length it encloses charge q. Therefore, from Gauss' theorem, 
iTirE = qlsQ 

or E = qjljis^r (1.6) 
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1.4 The differential form of Gauss' tlieorem 

Only a limited range of problems can be solved by the direct use of Equation (1.5). Another form, 
which is obtained by applying it to a small volume element, enables us to solve a much wider range of 
problems. Figure 1.4 shows such a volume element in Cartesian coordinates. 




Fig. 1.4 The elementary Gaussian surface used to derive the differential form of Gauss' theorem. 

To calculate the net flux out of the element, consider first the two shaded faces A and B which are 
perpendicular to the x-axis. The only component of E which contributes to the flux through these 
faces is E^. If the component of E is on A then, in general, the x component of E on 5 can be 
written + {dE^ / dx) dx . Provided that the dimensions of the element are small enough we can 

assume that these components are constant on the surfaces A and B. The flux of E out of the volume 
through the faces A and B is then 



E^ dy dz + 



dE 

E^ +^dx 
dx 



dE 

dy dz = — - dx dy dz 
dx 



The same argument can be used for the other two directions in space, with the result that the net flux 
of E out of the element is 



dE^ dE^ dE^ 



- + - 



- + - 



dx dy dz 



dx dy dz 



(1.7) 



Now if p is the local charge density, which may be assumed to be constant if the volume element is 
small enough, the charge enclosed in the volume is 



dQ = pdxdydz 



(1.8) 



17 



Electromagnetism for Electronic Engineers 



1. Electrostatics in free space 



Applying Gauss' theorem to the element and making use of Equations (1.7) and (1.8) gives the 
differential form of Gauss' theorem: 



dE^ dE dE^ 
+ — ^ + - 



dx dy dz 



— (1.9) 



The expression on the left-hand side of Equation (1.9) is known as the divergence of E. It is 
sometimes written as div E. The same expression can also be written as the dot product between the 
differential operator 



^.5 .5 . 5 ^ 
X hy hz — 

dx dy dz 



(1.10) 



and the vector E. In Equation (1.10) x, y , z are unit vectors along the x-, y- and z-axes. Using the 
symbol V , which is known as 'del'. Equation (1.9) can be written 

V-E = ^ (1.11) 

This abbreviation is not as pointless as it seems because Equation (1.1 1) is valid for all systems of 
coordinates in which the coordinate surfaces intersect at right angles. An appropriate form for V can 
be found for each such coordinate system. 

1.5 Electrostatic potential 

The electric field is inconvenient to work with because it is a vector; it would 
be much simpler to be able to work with scalar variables. The electrostatic 
potential difference (F) between two points in an electric field is defined as 
the work done when unit positive charge is moved from one point to the other. 



; 

Consider the figure on the right. The force on the charge is E, from Equation Charge Cu^^ 

(^ 'W Qn tViF» PYtprnnl fnrpf nPFrlprl tn VinlH it in F»niiili'Kri'iim i'q -TT. TViP» wnrV 



(1.3), SO the external force needed to hold it in equilibrium is -E. The work 
done on the charge by the external force when it is moved through a small 
distance dl is the product of the external force and the distance moved in the 
direction of that force. Thus the change in electrostatic potential is 



dV = -E A\ (1.12) 

The potential difference between two points A and B can be calculated by integrating Equation (1.12) 
along any path between them. Mathematically this is written 

Vb-V,=-\1^-A\ (1.13) 
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This kind of integral is called a line integral. This is a slightly tricky concept, but its application is 
limited in practice to cases where the symmetry of the problem makes its evaluation possible. 

The electrostatic potential is analogous to gravitational potential, which is defined as the work done in 
moving a unit mass against gravity from one point to another. The change in the gravitational 
potential depends only upon the relative heights of the starting and finishing points and not on the 
path which is taken between them. We can show that the same is true for the electrostatic potential. 
Figure 1.5 shows a possible path between two points A and B in the presence of the electric field due 
to a point charge Q at O. 




Fig. 1.5 When a charge is moved from >A to B in the field of another charge at O the change in 
electrostatic potential is found to be independent of the choice of the path APB. It depends only on the 
positions of the ends of the path. 
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The contribution to the integral of Equation (1.13) from a small movement dl of a unit charge at P is 
dV = dl 

4^^o^ (1.14) 

But the dot product r • dl is simply a way of writing 'the component of dl in the radial direction' 
using mathematical notation, and this quantity is just dr, the change in the distance from O. Thus 
Equation (1.14) can be integrated to give 



V -V =-[ — dr = 



(1.15) 



The potential difference between A and B therefore depends only on their positions and not on the 
path taken between them. By using the principle of superposition we can extend this proof to the field 
of any combination of charges. 

Potential differences are measured in volts. They are familiar to electronic engineers from their role in 
the operation of electronic circuits. It is important to remember that potentials are always relative. 
Any convenient point can be chosen as the zero of potential to which all other voltages are referred. 
Electronic engineers are inclined to speak loosely of the voltage at a point in a circuit when strictly 
they mean the voltage relative to the common rail. It is as well to keep this point in mind. 

It follows from the preceding discussion that the line integral of the electric field around a closed path 
is zero. This is really a formal way of saying that the principle of conservation of energy applies to the 
motion of charged particles in electric fields. In mathematical symbols the line integral around a 
closed path is indicated by adding a circle to the integral sign so that 

E-dl = 0 (1.16) 

The principle of conservation of energy often provides the best way of calculating the velocities of 
charged particles in electric fields. 

1.6 Calculation of potential in sinnple cases 

In simple cases where the electric field can be calculated by using Gauss' theorem it is possible to 
calculate the potential by using Equation (1.13). More complicated problems can be solved by using 
the principle of superposition. Since scalar quantities are much easier to add than vectors it is best to 
superimpose the potentials rather than the fields. 
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To show how this is done let us consider Fig 1.6 which shows a cross-sectional view of two long 
straight cylindrical rods each of radius a. The rods are parallel to each other with their centre lines d 
apart. Rod A carries a charge q per unit length uniformly distributed and rod B carries a similar charge 
-q. We wish to find an expression for the electrostatic potential at any point on the plane passing 
through the centre lines of the rods. 



Fig. 1.6 A cross-sectional view of a parallel-wire transmission line. 

The electric field of either rod on its own can be found by applying Gauss' theorem as described 
above with the result given by Equation (1.6). Since E is everywhere radial it follows that F depends 
only on r and 



V 



= -[ — - — dr = ^lnr + constant 



(1.17) 



where the constant of integration can be given any convenient value because the zero of potential is 
arbitrary. 

It is appropriate to choose the origin of coordinates to be at O, mid-way between the rods, because of 
the symmetry of the problem, so that OA lies along the x-axis. Then for rod A, r = \x - dl2 \ The 
same argument can be used for rod 5, giving an expression for the potential which is identical to 
(1.17) except that the sign is reversed and r = \x -\- d/2\. Superimposing these two results and 
substituting the appropriate expressions for the radii we get: 



V = - 



In 



2x -\- d 
2x - d 



(1.18) 



where the constant of integration has been set equal to zero. This choice makes V^O when x = 0. The 
same method could be used to find a general expression for the potential at any point in space. 
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1.7 Calculation of the electric field from the potential 

We have so far been concerned with means of calculating the potential from the electric field. In 
many cases it is necessary to reverse the process and calculate the field from a known potential 
distribution. Figure 1.7 shows how a small movement dl may be expressed in terms of its components 
as 

dl = xdx -\- y dy -\- z dz (1.19) 




Fig. 1.7 A small vector dl can be regarded as the sum of vectors dx, dy and dz along the coordinate 
directions 
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The electric field may likewise be expressed in terms of its components 

E = xE^+yE^+iE^ (1.20) 

Then from Equations (1.12), (1.19) and (1.20) the potential change along dl is given by 

dV = -{^E^+yEy+zE^)'{xdx + ydy + z dz) 

( . (1.21) 

= -\E^dx + Eydy +E^dz) 

Setting dy = dz = O^Q have, for a movement in the x-direction, 

dv 

Ex= 

dx 

with similar expressions for the other two coordinate directions. By superposition the total electric 
field is 



X hy hz 

dx dy dz 



(1.22) 



The expression in parentheses on the right-hand side of Equation (1.22) is termed the gradient of V. It 
can be obtained by operating on Kwith the operator V defined by Equation (1.10) so that Equation 
(1.22) may be written 

E = -gradF = -VF (1.23) 

This equation, like Equation (1.11), can be written in terms of any orthogonal coordinate system by 
using the appropriate form for W . 

From Equation (1 . 12) it can be seen that if dl lies in such a direction that V is constant it must be 
perpendicular to E. Surfaces on which Vis constant are known as equipotential surfaces or just 
equipotentials. They always intersect the lines of E at right angles. It has already been mentioned that 
field plots are useful aids to thought in electrostatics. They can be made even more useful by the 
addition of the equipotentials. Figure 1.8 shows, as an example, the field plot for the parallel wires of 
Fig. 1.6. 
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Fig. 1.8 The field pattern around a parallel-wire transnnission line 



1.8 Conducting nnaterials in electrostatic fields 



A conducting material in the present context is one which allows free movement of electric charge 
within it on a time scale which is short compared with that of the problem. Under this definition 
metals are always conductors but some other materials which are insulators on a short time scale may 
allow a redistribution of charge on a longer one. They may be regarded as conducting materials in 
electrostatic problems if we are prepared to wait for long enough for the charges to reach equilibrium. 
The charge distribution tends to equilibrium as exp(- t/r) , where the time constant i is known as the 
relaxation time. Some typical values are: 



Once the charges have reached equilibrium there can be no force acting on them and the electric field 
within the material must be zero. 

When an uncharged conducting body is placed in an electric field, the free charges within it must 
redistribute themselves to produce zero net field within the body. Consider, for example, a copper 
sphere placed in a uniform electric field. The copper has within it about 10^^ conduction electrons per 
cubic metre, and their charge is balanced by the equal and opposite charge of the ionic cores fixed in 
the crystal lattice. The available conduction charge is of the order of 10^^ C m"^, and only a tiny 
fraction of this charge has to be redistributed to cancel any practicable electric field. This 
redistribution gives rise to a surface charge, somewhat as shown in Fig. 1.9, whose field within the 
sphere is exactly equal and opposite to the field into which the sphere has been placed. 



copper 

distilled water 
fused quartz 



1.5 X 10"^^ s 
lOS 
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Fig. 1.9 The field pattern of the charge induced on a conducting sphere placed in a unifornn electric 
field. 
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This surface charge is known as induced charge. It is important to remember that the positive and 
negative charges balance so that the sphere still carries no net charge. The complete solution to the 
problem is obtained by superimposing the original uniform field on that shown in Fig. 1.9 to give the 
field shown in Fig. 1.10. Note that the flux lines must meet the surface of the sphere at right angles 
because the surface is an equipotential. 



Fig. 1.10 The field pattern around a conducting sphere placed in a uniform electric field. This pattern 
is obtained by superimposing the field of the induced charges (shown in Fig. 1.9) on the uniform field. 

Not only is there no electric field within a conducting body, but 
there is also no field within a closed conducting shell placed in 
an electric field. To prove this, consider the figure on the right, 
which shows a closed conducting shell Si. This must be an 
equipotential surface. If there is any electric field within Si there 
must be other equipotentials such as S2 lying wholly within Si. 
Now the interior of the shell contains no free charge so, applying 
Gauss' theorem to S2, the flux of E out of S2 is zero. But, since it 
has been postulated that S2 is an equipotential surface, this can be 
true only if E is zero everywhere on it and the potential of S2 is the 

A closed hollow earthed conductor can therefore be used to screen sensitive electronic equipment 
from electrostatic interference. The screening is perfect as long as there are no holes in the enclosure, 
for example to allow wires to pass through. Even when there are holes in the enclosure the screening 
can still be quite effective, for reasons which will be discussed in the next chapter. When the electric 
field varies with time other screening mechanisms come into play and the screening is no longer so 
perfect. 
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1.9 The method of images 

We have already seen that the electric field produced by a known distribution of charges can be 
calculated in simple cases, by the application of Gauss' theorem and the principle of superposition. In 
most practical problems, however, the charge is unknown and the problem is specified in terms of the 
potentials on electrodes. Simple problems of this type can be solved by the use of Gauss' theorem if it 
is possible to make assumptions about the distribution of charges from the symmetry of the problem. 

If an uncharged, isolated, conducting sheet is placed in an electric field, then equal positive and 
negative charges are induced on it. Normally this process requires currents to flow in the plane of the 
sheet, and the field pattern is changed so that the sheet becomes an equipotential surface. If, however, 
the sheet is arranged so that it coincides with an equipotential surface, the direction of current flow is 
normal to the plane of the sheet and the two surfaces become oppositely charged. If the sheet is thin, 
the separation of the positive and negative charges is small and the field pattern is not affected by the 
presence of the sheet. This fact can be used to extend the range of problems which can be solved by 
elementary methods. For example, a conducting sheet can be placed along the equipotential ^5 in Fig. 
1.8. It screens the two charged wires from each other so that either could be removed without 
affecting the field pattern on the other side of the sheet. Thus the field pattern between a charged wire 
and a conducting plane is just half of that of a pair of oppositely charged conducting wires. 
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The field between a charged wire and a conducting plane can be found by reversing the train of 
thought. We note that an image charge can be placed on the opposite side of the plane to produce a 
field which is the mirror image of the original field. The image charge is equal in magnitude to the 
original charge, but has the opposite sign. The plane is an equipotential surface in the field of the two 
charges, so it can be removed without altering the field pattern. The problem is then reduced to the 
superposition of the fields of the original and image charges. This method is known as the method of 
images. It can be applied to the solution of any problem involving charges and conducting surfaces if 
a set of image charges can be found such that the equipotentials in free space of the whole set of 
charges coincide with the conducting boundaries. 

1.10 Laplace's and Poisson's equations 

The method described in the previous section has been applied with ingenuity to a wide variety of 
problems whose solutions can be looked up when required. Unfortunately engineers are not free to 
choose the problems they wish to solve, and the great majority of practical problems cannot be solved 
by elementary methods. Figure 1.11 shows a typical problem: an electron gun of the kind used to 
generate the electron beam in a microwave tube for satellite communications. 



Cathode 



Electron 



Focus 
electrode 
at cathode 
potential ^ 




beam 



Equipotentials 



Fig. 1.11 The arrangement of a typical high-power electron gun. Such a gun might produce a 50 mA 
electron beam 2 mm in diameter for a potential difference between cathode and anode of 5 kV. 

In this case the field problem and the equations of motion of the electrons must be solved 
simultaneously because the space charge of the electrons affects the field solution. A general method 
which can be used, in principle, to solve any problem is obtained by combining Equations (1.9) and 
(1.22) to give 



d^V d^V d^V p 

— ^ + — 7 + — 

dx dy dzx 



(1.24) 
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This is known as Poisson's equation. It can also be written 

VV = -^ (1.25) 
where is given, in rectangular Cartesian coordinates, by 

V-^.^.il 0.26) 

ex' dy^ dz" 

When there is no free charge present the equation takes the simpler form known as Laplace's 
equation: 

c)^V d^V d^V 

W = ^ + ^ + ^ = 0 (1.27) 

dx' dy' dz' 

This equation has been solved for a very wide range of boundary conditions by analytical methods 
employing a variety of coordinate systems and by the special method known as conformal mapping, 
which applies to two-dimensional problems. These solutions can be looked up when they are required. 
Cases whose solutions are not available in the literature must, in nearly every case, be solved by 
numerical methods. When free charges are present in a problem it is necessary to use Poisson's 
equation as the basis of either an analytical or a numerical solution. There are only a few cases which 
can be solved analytically. 

In every kind of active electronic device electric fields are used to control the motion of charged 
particles. The methods described here can be applied to the motion of charged particles in vacuum. 
When the charge densities are small it is possible to calculate the electrostatic fields, neglecting the 
contributions of the charges to them, and then to integrate the equations of motion of the particles. At 
higher charge densities the fields are affected by the space charge and it is necessary to find mutually 
consistent solutions of Poisson's equation and the equations of motion. The motion of charge carriers 
in semiconductor devices such as transistors requires knowledge of the fields in material media as 
discussed in Chapter 2. 

1.11 The finite difference nnetliod 

The simplest numerical method for solving field problems is the finite difference method. In this 
method a regular rectangular mesh is superimposed upon the problem. The real continuous variation 
of potential with position is then approximated by the values of the potential at the intersections of the 
mesh lines. Figure 1.12 shows a small section of a two-dimensional mesh with a spacing h in each 
direction and the electrostatic potentials at the mesh points. 



29 



Electromagnetism for Electronic Engineers 



1. Electrostatics in free space 












— 1 


i 

( ^ 

1 

^ " 


A 






1 

L 

i 


1 






h 






* ^ » 





Fig. 1.12 Basis of the finite difference calculation of potential. 



To find an approximate relationship between the potentials shown we apply Gauss' theorem to the 
surface shown by the broken line. The component of the electric field normal to the section AB of the 
surface is given approximately by 

E„AB={K-V2)/h (1.28) 
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The flux of E through unit depth of the face AB is therefore 



= V -V 



(1.29) 



If the Gaussian surface does not enclose any charge the net flux of E out of it must be zero so that 



or 



4^0-^,-^2-^3-^4=0 



(1.30) 



Thus, if we know the potentials at points 1 to 4 approximately, we can use 
Equation (1.30) to obtain an estimate of Vq. Because the errors in the four 
potentials cancel each other out to some extent, and because the resulting error 
is divided by 4, the error in the value of Vq is normally less than the errors in 
the potentials used to calculate it. Equation (1.30) is conveniently summarized 
by the diagram on the right. 

This method can be used to find the fields around two-dimensional 
arrangements of electrodes on which the potentials are specified such as 
the concentric square tubes shown in the figure on the right. The method 
can be implemented on a spreadsheet as follows: 



a) A uniform square mesh is defined such that the electrodes coincide 

with mesh lines. The mesh spacing is chosen so that it is small enough to provide a reasonably 
detailed approximation to the fields whilst not being so small that the computational time is very 
large. 




b) Cells of the spreadsheet are marked out such that one cell corresponds to each mesh point. The 
symmetry of the problem can be used to reduce the number of cells required. Thus, for the 
geometry shown above it is sufficient to find the solution for one quadrant of the problem. 

c) The electrode potentials are entered into the cells corresponding to the electrodes and the formula 
in Equation (1.30) is entered into all the other cells. When symmetry has been used to reduce the 
size of the problem the formulae in the cells along symmetry boundaries make use of the fact that 
the potentials on either side of the boundary are equal. 

d) The formulae in the cells are then applied repeatedly ( a process known as iteration ) until the 
numbers in the cells cease to change. To do this the calculation options of the spreadsheet must be 
set to permit iteration. The final numbers in the cells are then approximations to the potentials at 
the corresponding points in space. 
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e) From this solution the equipotential curves can be plotted by interpolation between the potentials 
at the mesh points and the field components can be calculated at any mesh point. 

The method can be applied to more complicated problems including those with curved electrodes 
which do not fit the mesh and three-dimensional problems. Further information can be found in the 
literature. 

1.12 Summary 

In this chapter, starting from the inverse square law of force between two charges, we have derived a 
range of methods for solving practical problems involving electric fields in free space. The concepts 
of electric field, flux density and potential have been shown to be useful for these purposes. The ideas 
contained in this chapter find their direct application in problems about voltage breakdown between 
electrodes in air and those dealing with the motion of charged particles in vacuum. 

The very limited range of problems which can be solved by elementary methods can be extended by 
the use of the principle of superposition and the method of images. In most real problems, however, 
the electric field can be calculated only by solving Laplace's or Poisson's equations. Cases which 
have not been solved before generally have to be tackled using numerical methods. 
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2. Dielectric materials and capacitance 

Objectives 

• To discuss how and why an electric field is affected by the presence of dielectric materials 

• To introduce the electric flux density vector D as a useful tool for solving problems involving 
dielectric materials. 

• To derive the boundary conditions which apply at the interface between different dielectric 
materials. 

• To introduce the idea of capacitance as a general phenomenon which is not restricted to 
capacitors. 

• To demonstrate the calculation of capacitance by the use of Gauss' theorem, field solutions 
and energy methods. 

• To introduce the idea of stored energy density in an electric field. 

• To discuss the causes of electrostatic interference and techniques for reducing it. 
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2.1 Insulating nnaterials in electric fields 



Very many materials do not allow electric charges to move freely within them, or allow such motion 
to occur only very slowly. These materials are not only used to block the flow of electric current but 
also to form the insulating layer between the electrodes of capacitors. In this context they are known 
as dielectric materials. By making an appropriate choice of dielectric material for a capacitor it is 
possible to reduce the size of a capacitor of given capacitance or to increase its working voltage. 

If a dielectric material is subjected to a high enough electric field it becomes a conductor of 
electricity, undergoing what is known as dielectric breakdown. This controls the maximum working 
voltage of capacitors, the maximum power which can be handled by coaxial cables in high-power 
applications such as radio transmitters, and the maximum voltages which can be sustained in 
microcircuits. It is not always appreciated that because dielectric breakdown depends on the electric 
field strength it can occur when low voltages are applied across very thin pieces of dielectric material. 

In order to understand the behaviour of dielectric materials in electric fields it is helpful to make a 
comparison with that of conductors. Figures 2.1 (a) and (b) show respectively a conducting sheet and a 
dielectric sheet placed between parallel electrodes to which a potential difference has been applied. 
The potential difference is associated with equal and opposite charges on the two electrodes. 
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(a) ib) 
Fig. 2.1 Comparison between (a) a conductor and (b) an insulator placed in an electric field. 



The conducting sheet of Fig. 2.1 (a) contains electrons which are free to move and set up a surface 
charge which exactly cancels the electric field within the conductor in the manner discussed in 
Chapter 1 . The electrons in the dielectric material, on the other hand, are bound to their parent atoms 
and can only be displaced to a limited extent by the applied electric field. This displacement, however, 
is sufficient to produce some surface charge and the dielectric is then said to be polarized. The 
surface charge is not sufficient to cancel the electric field within the sheet, but it does reduce it to 
some extent, as shown in Fig. 2.1 (b). 



Polarization may also produce a volume distribution of charge, but we shall assume that this does not 
occur in the materials in which we are interested. It is important to remember that the surface charge 
produced by the polarization of a dielectric is a bound charge which, unlike the surface charge 
induced on a conductor, cannot be removed. The polarization charges must also be carefully 
distinguished from any free charge which may reside on the surface of a dielectric. 
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On materials which are good insulators, free charges may persist for long periods and strong electric 
fields may build up as a result of them. These phenomena have many important practical 
consequences, but they are not easy to study theoretically because the distribution of the charges is 
usually unknown. In metal oxide semiconductor (MOS) integrated circuits, for example, it is possible 
for charges to build up on the gate electrodes if they are left unconnected. The electric field produced 
by these charges can be strong enough to cause dielectric breakdown of the silicon dioxide layer. This 
is why special precautions have to be taken when handling these circuits. In what follows we shall 
assume that the dielectric is initially uncharged and that any surface charge is the result of 
polarization. 

To put this subject on a quantitative basis, let us suppose that the electrodes in Fig. 2.1 (b) carry a 
surface charge a per unit area and that the surface charge on the dielectric is Op per unit area. Now, 
assuming that the electric field is everywhere uniform and normal to the electrodes, it can be shown 
that the field outside the dielectric is given by 

Ea=- (2.1) 

This result is obtained by applying Gauss' Theorem to the field between plane, parallel, electrodes. 
Similarly the field within the dielectric is 

E,= (2.2) 

Equation (2.2) can be rewritten as 

(2.3) 



^0 ^O^r 



where = aji^ — a is known as the relative permittivity of the material. Since cr^ < a it 

follows that Sr > 1. It is unfortunate that the symbol Sr has been adopted for this property of dielectric 
materials because there could be some confusion between it (a dimensionless quantity) and the 
permittivity, defined as e = soSr, and measured in Farads per metre. Care must be taken not to get 
these symbols confused with each other. 

In order to make the theory simpler, we shall assume that ap, is proportional to a and that Sr is 
therefore a constant. This assumption holds good for many of the materials used in electronic 
engineering, but it is very important to remember that it is not always valid. In particular, for some 
materials, Sr may depend on: 



35 



Electromagnetism for Electronic Engineers 



2. Dielectric materials and capacitance 



• the strength of the electric field; 

• frequency (if the field is varying with time); 

• the orientation of crystal axes to the field; 

• the previous history of the material. 

Problems involving linear dielectric materials could be solved by calculating the polarization charges 
and finding the fields resulting from both the free charges on the electrodes and the bound 
polarization charges. This would not usually be easy and it is much better to use an approach in which 
the polarization charges are implicit. To do this we introduce a new vector known as the electric flux 
density, which is defined by 



In the example given above the electric flux density outside the dielectric isD = a and that within the 
dielectric is likewise a. In other words D depends only on the free charges, unlike E, which depends 
on the polarization charges as well. 

It can be shown that, subject to the validity of the assumption that £r is a constant, the argument given 
above can be generalized to cover pieces of dielectric material of any shape. Gauss' theorem 
(Equation (1.5)) can thus be written in a form which is valid for problems involving dielectric 
materials: 



D = ^E 



(2.4) 
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(2.5) 



or, in differential form 



V B = p 



(2.6) 



2.2 Solution of problems involving dielectric nnaterials 

Many problems in electrostatics deal with sets of electrodes together with dielectric materials. When 
the symmetry of a problem is simple it is possible to use Gauss' theorem in much the same way as in 
Chapter 1. 

Figure 2.2 shows a coaxial cable in which the space between the conductors is filled with a dielectric 
material of permittivity s. We wish to find an expression for the electric field within the dielectric 
when the potential difference between the electrodes is Vq. 




Fig. 2.2 The arrangement of coaxial cylindrical electrodes, an idealization of a coaxial cable. 

Assume that the conductors carry charges ± q per unit length, with the inner conductor being 
positively charged. Applying Gauss' theorem as given in Equation (2.5) to a cylindrical Gaussian 
surface of radius r we find that the radial component of D is given by 



27rr 



(2.7) 



then, using Equation (2.4) 




27rsr 



(2.8) 



The potential difference between the conductors is found by using Equation (1.13) 




(2.9) 
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Then 



Ins 



(2.10) 



The charge q which was assumed for the purposes of the solution need not be calculated explicitly and 
can be eliminated to give the radial component of the electric field as a function of radial position: 

V 

E = ^ r (2.11) 

r\n\bla) 

2.3 Boundary conditions 

When two or more dielectric materials are present it is necessary to treat each region separately and 
then to apply the appropriate boundary conditions at the interfaces. There are three of these 
conditions relating to V, E and D respectively. 

The electrostatic potential is continuous at a boundary, that is, it cannot change suddenly there. The 
physical reason for this condition is that an abrupt change in the potential would imply the presence of 
an infinitely strong electric field. 



To find the boundary condition for the electric field we 

consider an infinitesimal closed path as shown in the figure on *.';>.* '' V* • • 

the right. The path is chosen so that it crosses the boundary 
between two dielectrics having permittivities and S2 as 
shown. If the loop is made very thin, then the contributions to 
the line integral in Equation (1.16) arising from the parts of the 
loop normal to the boundary are negligible. If, in addition, the 
tangential components of the electric field are Eti and Et2, then the integral becomes 

{E,,-E,,)dl = 0 
SO that 

En=E,2 (2.12) 
This result can be stated in words as: the tangential component ofK is continuous at a boundary. 



The boundary condition for the electric flux density can be 
found in a similar way by using Gauss' theorem. The figure in 

the margin shows a boundary between two dielectric materials ] : ^ \ 1 Dn 



with an infinitesimal Gaussian surface which crosses it. If the 
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thickness of the 'pill box' is very small, then the contributions to the flux from the parts of the surface 
which are normal to the boundary are negligible and the integral becomes 

{D„,-DjdA = 0 

where Z)„i, and Z)„2 are the components of D normal to the boundary on each side of it and dA is the 
area of the part of the boundary lying within the Gaussian surface. Thus 

D„,=D„, (2.13) 
or, in words, the normal component ofD is continuous at a boundary. 

To solve problems with two, or more, layers of dielectric material we first apply Gauss' theorem to 
find D everywhere since that does not depend on the presence of the materials. Next Equation (2.4) is 
used to find E in each region. Finally the potential difference across each layer is found using 
Equation (1.13). 
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2.4 Capacitance 



Capacitors are very familiar as circuit elements, but it is not always realized that the idea of 
capacitance is more general. Capacitance exists between any pair of conductors which are electrically 
insulated from each other. Thus there is a capacitance between adjacent tracks on a printed circuit 
board, but this does not usually appear in the circuit diagram. This 'stray' or 'parasitic' capacitance 
can cause unwanted coupling between the parts of a circuit, causing it to oscillate or misbehave in 
some other way. Very few electronic engineers ever need to calculate the capacitance of a capacitor; 
they are much more likely to need to estimate the magnitude of a stray capacitance. 

The figure on the right shows a cross-sectional view of two 

adjacent tracks on a printed circuit board. Let us suppose that the ^. q _ q 

A f ^ I 

tracks are insulated from each other and from earth and that < y / / / / 



charge Q is transferred from B\o A. Electrode B must then carry L. Z L. €- 1. — i 

charge -Q. As a result of this transfer of charge, an electric field 

exists around the electrodes such that the potential difference between them is V. 

If the dielectric material of the printed circuit board has a permittivity which does not vary with the 
electric field, then the system is linear and the principle of super-position may be applied. It follows 
that the potential difference between the electrodes is directly proportional to the charges on them, so 
we can write 



where C is a constant of proportionality which is readily recognized as the capacitance between the 
electrodes familiar from elementary circuit theory. The unit of capacitance is the farad (F) and 



Stray capacitances are usually of the order of picofarads. 

2.5 Electrostatic screening 

It has already been noted that unwanted capacitive coupling between electronic circuits can be a major 
problem. This is part of the larger problem of electromagnetic interference; another aspect, inductive 
coupling, is discussed in Chapter 6. The problem with all types of electromagnetic interference is how 
to minimize it rather than how to calculate its magnitude accurately. Electromagnetic theory provides 
the means for understanding the causes of electromagnetic interference and the techniques for dealing 
with them. 

A simple case of the coupling of two circuits by stray capacitance is shown in Fig. 2.3. The circuits 1 
and 2 are linked by the stray capacitance C and by a common earth. The stray capacitance is small, 
typically of the order of a pico-farad, so its impedance (Z^ =1/ jcoC^ ) is high, but decreases with 
increasing frequency. In this problem V\ is the source of the interference picked up by circuit 2. 




Q = CV 



(2.14) 



1 F = 1 C V" . Most capacitances are small and measured in microfarads, nanofarads or picofarads. 
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Fig. 2.3 Showing how two circuits can be coupled together by a stray capacitance Cg. The unwanted 
signal coupled fronn circuit 1 into circuit 2 can be large enough to be troublesonne. 

The current flowing through the capacitor is small compared with that in so the spurious signal 
appearing at the input of the amplifier is approximately 



R 



Ll 



(2.15) 



where R'.^ = {Rs2^m )/ (^^2 + ^in ) • The spurious voltage given by Equation (2. 1 5) is to be compared 
with the signal voltage at the input of the amplifier 

K.,=V2Rj{Rs2+R,n) 



Equation (2.15) shows that the spurious signal is greatest when the source impedance R^u is low and 
the effective input impedance of the second circuit R'in is high. An indication of the order of 
magnitude of capacitance which can cause trouble can be obtained by supposing that circuit 1 is the 
a.c. main. R^i is then very small. If R\n is 1 MQ and is 1 |iV, Zs is approximately 2.4 x 10^^ Q, 
which corresponds at 50 Hz to a stray capacitance of the order of 10"^"^ F. 

The capacitive coupling between the circuits can be reduced by putting an earthed screen between 
them, as shown in Fig. 2.4. The stray capacitance is divided into two parts in series with each other 
and with their common point earthed. In practice, unless the screen completely encloses one of the 
circuits, there is still a residual capacitance connected directly between P and Q, bypassing the screen. 
Quite effective screening can be achieved with a partial enclosure provided that it intercepts most of 
the field lines passing from P to Q. It is essential that the screen is earthed, otherwise P is connected 
to Q by Csi in series with Cs2, without a signal path to earth from their common point. 
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Fig. 2.4 The unwanted coupling between the two circuits illustrated in Fig. 2.3 can be greatly reduced 
by putting an earthed screen between thenn. Any current passing through Csi is conducted to earth 
instead of passing through Cs2 to O. 

When very good screening against low-frequency electric fields is required then a closed conducting 
box must be used. Since all the sources of the field lie outside the box and it must be an equipotential 
surface it follows that, theoretically, the electric field inside it is zero. Such a box is sometimes known 
as a 'Faraday cage'. In practice, any enclosure will have one or more holes in it to allow wires to pass 
in and out and these may reduce the effectiveness of the screen. Similarly, the screening effectiveness 
is affected by the way in which the joints of the box are made. 
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The electric screening effectiveness of an enclosure is defined as the ratio of the magnitude of the 
electric field with the screen to that at the same point when the screen is removed. It is usually 
expressed in decibels. At high frequencies other factors come into play. 

The effects of capacitive coupling can also be reduced by using a differential amplifier if the source 
V2 can be isolated from earth. Figure 2.5 shows how this works. The stray capacitances Csi and Cs2 are 
often approximately equal and so the spurious signals appearing at the normal and inverting inputs of 
the amplifier are also very nearly the same. If the amplifier has a high common-mode rejection ratio, 
then only their difference is amplified and added to the wanted signal. Adding an earthed screen 
reduces the stray capacitive coupling and therefore reduces the unwanted signal still further. The 
source for circuit 2 may be a transducer situated some distance from the amplifier. In that case the 
connecting cables must be screened as well as the amplifier, a point which is discussed further in 
Chapter 6. 




Fig. 2.5 When the input to a circuit can be isolated from earth it is possible to reduce capacitive 
coupling effects by using an amplifier with a differential input. The unwanted signals coupled through 
Csi and Cs2 are approximately equal but only their difference is added to the signal. 

It is also possible for there to be capacitive coupling between the output and the input of an amplifier. 
When that happens there can be positive feedback causing the circuit to oscillate, so once again 
screening is needed. 

2.6 Calculation of capacitance 

In simple cases where the field problem can be solved by the direct application of Gauss' theorem it is 
straightforward to calculate the capacitance. The sequence of steps is: 

1 . Assume charges ± ^ on the conductors. 

2. Apply Gauss' theorem, superposition or the method of images to find E. 

3. Integrate E along any convenient path between the conductors to find the potential difference 
between them. 

4. Calculate the capacitance from C = q /V. 

When the shapes of the electrodes are more complicated than those in the examples above it is no 
longer possible to use the same method to calculate the capacitance. An alternative method is: 
1. Assume potentials 0 and Von the electrodes. 
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2. Solve Laplace's equation analytically or numerically to give values of the potential 
everywhere. 

3. Calculate E close to the surface of one of the electrodes using E = -grad V. 

4. Calculate the surface charge density distribution on the electrode using a = sE. 

5. Integrate the charge density distribution over the surface of the electrode to give the total 
charge q. 

6. Calculate the capacitance from C = q / V. 

To summarize: capacitance can be obtained by calculating either 

a) The potential difference between the electrodes from assumed charges, or 

b) The charges on the electrodes from an assumed potential difference between them. 

By using both approaches in an approximate way it is possible to obtain upper and lower bounds to 
the capacitance. The formal application of this in the energy method is described later in the chapter. 

2.7 Energy storage in the electric field 

When charge is transferred from one plate of a capacitor to the other work is done against the electric 
field. This work is stored as potential energy in the capacitor. It is easy to show that the stored energy 
is given by 



It is important to remember that Equation (2.16) applies to any pair of electrodes with a potential 
difference between them, not just to capacitors as lumped components. 



Fig. 2.6 The electric field between two charged conductors can be thought of as being made up of 
tubes of flux like the one shown. The walls of the tube are everywhere parallel to the local direction of 
the electric field. Its ends are terminated by equal and opposite charges. 

It is sometimes useful to think of the energy stored in a capacitor as being distributed throughout the 
electric field associated with it. Figure 2.6 shows a capacitor with electrodes of arbitrary shape. We 
define a flux tube by considering a small element of the surface of the conductor A and following the 
lines of E which start from its boundary through space until they terminate upon B. Since the flux 
lines can never meet or cross each other it follows that the whole electric field of the capacitor can be 
divided up into flux tubes. In addition we observe that, if the charge on the element of A from which 



W = -CV' = 
2 



2 C 2 




(2.16) 
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the tube starts is +dQ, then that on the element of B on which it ends must be -dQ by a straightforward 
apphcation of Gauss' theorem to a Gaussian surface enclosing the tube. 

The figure on the right shows a small part of the flux tube on a larger 
scale. A short length of the tube is chosen having length dl as shown, 
and the cross-sectional area of the tube at this point is dA. If the local 
electric field strength is E, then the potential difference between the 
ends of this element of volume isdV = E dl, where the dot product can 
be omitted because E and dl are parallel to each other by definition. 
Likewise, the relationship between the electric flux in the tube and the 
charges at its ends means that dQ = D dA. 

The energy stored in the element of volume is then 

dW = -dOdV = -DEdldA = -DEdv 
2 2 2 

using Equation (2.16) with the volume of the element represented by dv. In the limit as dv shrinks to 
zero the energy density in the electric field is /^DE. A more rigorous derivation of this result which 
allowed for the possibility that the vectors D and E are not parallel to each other would show that the 
energy density in the field is given by 
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w = -D-E (2.17) 
2 

and the total energy stored in the field of the electrodes is 

fF^ijjj D-E dv (2.18) 
where the integral is taken over all that part of space in which the electric field is not zero. 

2.8 Calculation of capacitance by energy nnethods 

In many cases the solution required to a field problem is not detailed field information but a single 
number, the capacitance. The use of the concept of stored energy provides a particularly simple way 
to obtain estimates of capacitance without finding the exact solution to the field problem. The starting 
point is the fact that when a physical system is in stable equilibrium its stored energy is normally a 
minimum. From this it follows that any perturbation of the system must result in an increase in the 
stored energy. 

Consider now two possible perturbations of the system of charges and electric fields in a capacitor. In 
the first case we retain the electrodes as equipotential surfaces but alter the charges on them in such a 
way that the equipotentials in the inter-electrode space assume particularly simple shapes. From these 
perturbed equipotentials we calculate the energy density in the field and, hence, the stored energy W 
of the system. This energy must be greater than the stored energy W of the unperturbed system, so 
that 

1 2W' 
W'>-CV^ or (2.19) 
2 

This provides an upper bound to the capacitance. 

In the second case we hold the charges on the electrodes constant but redistribute them in such a way 
that the flux lines have a particularly simple form. When this is done the electrodes are no longer 
equipotential surfaces. Once again the stored energy, say W", is calculated, and this time 

W">— or C>^— (2.20) 
2C 2W" 



giving a lower bound for the capacitance. 
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2.9 Finite element nnethod 

Many computer packages for the solution of field problems are based upon the use of finite elements. 
This technique is not suitable for hand calculation and the algorithm is a little complicated so only the 
principles of the method as applied to two-dimensional problems are described here. 

In the finite element method the region between the electrodes is divided into a large number of 
triangular elements. Figure 2.7(a) shows one such element in the x-y plane. It can be shown that the 
condition that Laplace's equation should be obeyed inside the element is equivalent to the 
requirement that the energy stored within it should be a minimum. If the potential within the element 
is approximated by 

V{x,y) = a-\-bx-\-cy (2.21) 
then the electric field components are 

E^=-b and E^=-c (2.22) 
and the energy density is 

^ = ^^^{t^ +c^) (2.23) 

Now the coefficients b and c can be expressed in terms of the positions of the vertices of the triangle 
and of the potentials at them. Thus the stored energy in the element can be expressed in terms of the 
same quantities. 





fa) f (b) 

Fig. 2.7 The finite element method: (a) a single element and (b) a cluster of elements. 



Now consider the cluster of triangular elements with a common vertex at^, as shown in Fig. 2.12(b). 
The energy stored in the cluster of elements, W, is the sum of the energies stored in each. We now 
select the value of the potential at A which minimises the stored energy by requiring that 
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dW 



= 0 



(2.24) 



The result is an equation giving the potential at A in terms of the potentials at the other vertices and 
their positions. It is interesting that, if triangular elements are formed from the square elements in Fig. 
1.12 by adding diagonals from the lower left to the upper right comer of each square, then the 
application of the finite element method yields Equation (1.30). 

The finite element method is important chiefly because of the freedom with which the sizes and 
shapes of the elements may be chosen. It is, therefore, possible to select a set of elements which fits 
the boundaries of the problem with good accuracy and which provides a concentration of small 
elements in regions where the potential is changing rapidly. In two-dimensional problems the spatial 
variation of the potential can be represented as a smooth surface. The finite element method 
approximates this surface by one with many small flat facets rather like the surface of a cut gem. 
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2.10 Boundary element nnethod 

The finite difference and finite element methods both depend upon the division of the space between 
the electrodes into a large number of small elements. The potential values at the nodes of the mesh are 
manipulated according to some algorithm until the solution has converged. The result of the process, 
in each case, is detailed information about the field at every point within the solution space. In large 
three-dimensional problems the number of variables can be very large indeed and the computation 
time correspondingly long. 

Very often the result required from a field calculation is not the detailed field information but only a 
single number such as the capacitance. For problems of this kind an alternative method, the boundary 
element method, is useful because it only requires mesh points and potentials to be specified on the 
boundaries. This has the benefit that the number of variables in the solution is less than in the 
corresponding finite difference or finite element solution by two or three orders of magnitude. It also 
makes it easier to model problems with complicated boundaries. The fields at internal points can be 
calculated from the solution at the boundary if required. 

2.11 Sunnnnary 

Insulating materials are widely used in electronic engineering both to provide electrical isolation and 
to enhance the performance of capacitors. The dielectric qualities of most insulating materials can be 
represented by a constant, the permittivity. It is convenient to introduce a new field vector, the electric 
flux density D, whose value is independent of the presence of dielectric materials. The use of this 
vector makes the solution of problems involving dielectric materials easier. It is important to use the 
correct boundary conditions at the interfaces between different dielectric materials. 

Capacitance occurs in capacitors and in stray capacitance between conductors. It can be calculated 
either by a direct solution of the field problem or by considering the stored energy in an 
approximation to the field associated with the capacitor. The latter method gives quite good accuracy 
with relative ease. Stray capacitances lead to unwanted coupling between electronic circuits. 
Techniques for reducing the coupling can be understood using the theory of electrostatics. 
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3. Steady electric currents 



Objectives 

• To introduce the ideas of current flow in electrical conductors and of electrical resistance. 

• To show that Ohm's law is a special case of a more general relationship between current 
density and electric field. 

• To derive expressions for the power dissipated per unit volume within a conductor. 

• To show how the distribution of current within a conductor can be calculated. 

• To introduce the continuity equation and to show how Kirchhoff s current law is a special 
case of it. 

• To introduce the concept of the electromotive force of a source of electric power. 

• To modify the equation for the line integral of the electric field around a closed path derived 
in Chapter 1 to allow for the presence of sources of electric power, and to show that the 
resulting equation is a generalization of Kirchhoff s voltage law. 

• To show how the resistance of a conductor can be calculated directly, by the use of Laplace's 
equation, and by energy methods. 



In the preceding chapters we have studied problems involving charges which are at rest. This is 
convenient for the formal exposition of electromagnetic theory, but most people meet electrical 
phenomena for the first time in the form of electric currents flowing in wires. A current can flow in 
any medium in which there are charges which are free to move. These conduction charges may be 
electrons, positively charged 'holes', or positive or negative ions, according to the material. They are 
in continuous random motion, with a distribution of velocities which depends upon the temperature of 
the material. They are also constantly colliding with each other and with the atomic structure of the 
material. 

When a conducting material is placed in an electric field the conduction charges are accelerated in the 
direction of the field. The velocity acquired is small compared with the average value of the random 
velocity at ordinary temperatures, and is superimposed on it. The ordered part of the motion would 
increase without limit were it not for the collisions which convert it into random motion. Overall, the 
effect of the field is to add a small average drift velocity to the random velocity. The magnitude of the 
drift velocity is related to the strength of the field by the equation 



where fi is known as the mobility of the charge carriers. The mobility depends on 

1 . the type of charge carrier; 

2. the material in which it moves; 

3. the temperature of the material. 
It can also depend on 



3.1 Conduction of electricity 



(3.1) 
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4. the strength of the electric field; 

5. the orientation of the crystal axes to the field. 

For many materials // is a scalar quantity which varies only with temperature and this will be assumed 
here. 

The rate of flow of electric charge across unit area of a plane normal to the direction of is the 
current density given by 



where n is the density of charge carriers and q the charge on each. Combining Equations (3.1) and 
(3.2) gives 



where a is the conductivity of the material, measured in Siemens per metre. J is measured in 
coulombs per square metre per second or in amperes per square metre, defining the ampere as a 
current of one coulomb per second. It is sometimes convenient to use the reciprocal of conductivity 
which is known as resistivity (p). Resistivity is measured in V m A"^ or Q m. 



J = nq\^ 



(3.2) 



J = nqjuE = (jE or E = pJ 



(3.3) 
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To link these ideas with those of elementary circuit theory 
consider the flow of electric current in a section of straight 
wire as shown in the figure on the right. If the wire has cross- 
sectional area A and carries a current / uniformly distributed 
across it, then the current density is 



^ Cross- sectional 
area A 



d 



J = I/A 



(3.4) 



The relationship between the electric field in the wire, the potential difference between its ends, and 
its length is 



E = V/l 



(3.5) 



But E = pj from Equation (3.3) and therefore 



V = 



A 



I = RI 



(3.6) 



which is the familiar form of Ohm's law. Equation (3.3) is thus shown to be a general form of Ohm's 
law applicable to problems in which the current density is not uniform. 

3.2 Ohmic heating 

When conduction charges are accelerated by an electric field, they continually gain energy from it. 
This ordered kinetic energy is transferred to the bulk of the material by collisions, so increasing the 
random thermal motions of the atoms. 

The conversion of electrical energy into heat is familiar from its everyday use in electric heaters and 
light bulbs. For a steady current / sustained by an applied voltage Fthe power input is 



p = vi = i^r = v^Ir 



(3.7) 



an equation which is well known from elementary circuit theory. By using the results of the previous 
section we can find corresponding formulae for the details of the energy dissipation within a 
conducting material. 



Consider the element of volume shown in the figure on the right. Provided 
that the properties of the material are the same in all directions, we can set 
the X-axis parallel to the direction of the current without loss of generality. 
The potential difference across the element is 
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V = ESx (3.8) 
and the current flowing though it is 

I = JSy& (3.9) 
Thus the power dissipated in the element is 

SP = EJSxSy& = EJSv = aE'Sv = pj'Sv (3.10) 

where Sv is the volume of the element. A more rigorous derivation would show that 

SP = E'3Sv (3.11) 

This expression is valid even when the vectors E and J are not parallel to each other. Equation (3.10) 
shows that the power density is greatest in regions of high current density. This phenomenon is of 
great practical importance. In fuses the provision of a section of thin wire in a circuit ensures that that 
section is heated to melting point before any other part of the circuit is damaged. On the other hand, 
wiring joints which are carelessly designed or made can overheat and fail. 

3.3 The distribution of current density in conductors 

We would like to know how to calculate the variation of current 
density with position. To do this we consider the net current 
flow out of a closed surface S enclosing a volume Fas shown in 
the figure on the right. If the current density in the region of S is 
J (not necessarily a constant) then the net current flow out of S 
is given by the integral of the normal component of J over S, 
that is 

current out of 5 = J • dA (3.12) 



Now if there is a net current flow out of a closed surface, then the total charge enclosed by that 
surface must be changing with time because charges cannot be created or destroyed. Thus 




d 

current out of S = (charge enclosed hy S) 

dt 



which may be written in mathematical notation as 
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§J-dA = -|jnP^v (3.13) 



This is a mathematical statement of the principle of conservation of charge. Note that p is charge 
density, not resistivity, in this equation. 

Using the same method as we used in Chapter 1 to get from Equation (1.5) to Equation (1.11), it is 
easy to show that the differential form of Equation (3.13) is 

V-i = (3.14) 

ot 

This equation, in either its integral or its differential form, is known as the continuity equation. It is 
valid for the flow of charges in a vacuum as well as in conducting materials. 

When the current flow is in a steady state, the right-hand sides of Equations (3.13) and (3.14) are zero. 
If, furthermore, the conductivity of the material is constant, then 

V-J = V-(crE)=cr(V-E)=0 (3.15) 
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since cr is a constant. But 



E = -VV 



(Equation (1.23)) 



so that, from Equation (3.15) 

Thus when a steady electric current flows in a material of constant conductivity the potential 
distribution obeys Laplace's equation (Equation 1.27). 

In the special case of currents flowing in wires, as shown in Fig. 3.1, the integral on the left-hand side 
of Equation (3.13) becomes the sum of the currents in the wires 



(3.16) 



which is Kirchhoff's current law familiar from electrical circuit theory. The continuity equation is 
therefore a generalization of Kirchhoff s current law. 




Fig.3.1 Applying the continuity equation to the special case of steady currents flowing in wires 



In general, provided that the conductivity is constant, the current 
density distribution can be calculated from the solution of Laplace's 
equation. The figure on the right shows a typical problem: a strip with 
a step change in its width. In this case Laplace's equation would have 
to be solved by numerical methods. The potential is specified by 
assuming that at some distance from the step the equipotentials are 

straight lines perpendicular to the length of the strip. Along the sides there is a boundary between a 
conductor and an insulator. Since no current can cross this boundary the electric field must be 
tangential to the boundary. 
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3.4 Electric fields in the presence of currents 

When electric currents flow through a system of conductors there is an electric field in the space 
around the conductors in addition to the field within the conductors which drives the current. As an 
example, consider the field around a coaxial line carrying equal and opposite currents in the core and 
the sheath, as shown in Fig. 3.2. Within the conductors the field lines, like the current fiow lines, lie 
parallel to the axis. Outside them the field is modified fi-om the purely radial electrostatic solution. 

The finite resistance of the conductors requires that the potential at A should be higher than that at B 
and that at C should, likewise, be greater than that at D. There must, therefore, be axial components of 
the electric field in the space between the conductors as well as in the conductors themselves. The 
tangential component of the electric field is continuous at the boundary (Equation (2.12)) so the field 
lines are curved as shown and are no longer normal to the boundary. 



/ D 




c 














1 

B 















///////////// 



Fig. 3.2 Electric field lines and equipotential surfaces in a coaxial cable with resistive conductors. 



3.5 Electronnotive force 



In this chapter we have assumed, so far, that sources of electric current exist, but we have paid no 
further attention to them. These sources have been represented as idealized voltage generators. In 
practice, they are most likely to be electrochemical cells or electromagnetic generators. The details of 
the inner workings of these sources do not fall within the scope of this book, but it is important to 
consider the effects of sources in the theory of electromagnetism. 

We shall regard a source as a device capable of maintaining a potential difference and a steady fiow 
of electric current between its terminals. Such a device is a source of electric power because it does 
work in driving the electric current around the circuit against the circuit resistance. As the charges 
move around the circuit from one terminal of the source to the other, they lose potential energy as it is 
converted into heat. 



In Chapter 1 we saw that in an electrostatic field the line integral of the electric field around a closed 
path is zero 
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This equation remains true near a current-carrying circuit provided that the path of integration does 
not pass through the generator. But if the line integral of the electric field is evaluated around the 
circuit from one terminal of the source to the other, the result is not zero. This is because a net amount 
of work has to be done to move the charges along this path against the resistance of the circuit. The 
potential energy given up by a charge in this way is restored when it moves through the source. This 
can be expressed mathematically by writing 



where is the electromotive force (e.m.f.) of the source and the integral is taken around the circuit 
(Q from one terminal of the source to the other. 

When the source is imperfect it is usual to equate the electromotive force with the Thevenin voltage, 
i.e. the open-circuit voltage of the source, and to include the source impedance in the external circuit. 




(3.17) 



c 
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For the special case of a circuit consisting of lumped components and generators joined by lossless 
conductors, the integral in Equation (3.17) becomes the sum of the potential differences across the 
components. The electromotive force is, likewise, the sum of the e.m.f 's of all the generators in the 
circuit. The result is Kirchhoff's voltage law: 



4- R T = 4- 4- 



(3.18) 



3.6 Calculation of resistance 



From Ohm's law in the form 



V = 



Pl 
A 



I = RI 



(Equation (3.6)) 



it is clear that resistances can be calculated either by assuming a current and working out the potential 
difference, or by the opposite process. It is also possible to make use of energy methods, with the upper 
and lower bounds for the resistance being given by these two approaches. In simple cases the resistance 
of a conductor can be calculated directly using the expression in the brackets in Equation (3.6). 

This approach is only of use in the relatively small number of cases where the shape of the conductor 
is such that either the current, or the potential, distribution is simple. This is clearly not the case with 
shapes such as that shown in Fig. 3.3. Problems of this kind can be solved by the following procedure: 

1. Assume a value for the potential difference F between the ends of the strip; 

2. Solve Laplace's equation with the appropriate boundary conditions; 

3. Calculate the electric field and current density distributions on any convenient line across the 
strip using E = -V V and / = aE ; 

4. Integrate J across the strip to find the total current /; 

5. R = V/I. 

These steps can be carried out by numerical methods if necessary. 
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Fig. 3.3 Electric field lines (shown dashed) and equipotential surfaces for a current flowing through a 
conducting strip with a step change in width. 



This method can sometimes be a bit tricky to apply because of the need to ensure that the normal 
component of E is zero at the edge of the strip. We can sometimes get round the difficulty by using 
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the principle of duality. In a field map such as Fig. 3.3 the field lines and the equipotentials are 
always at right angles to each other. We can therefore imagine a second problem whose field pattern 
is identical to the first except that the roles of the field lines and the equipotentials have been 
exchanged. Two problems related to each other in this way are called duals of each other. Sometimes 
it is easier to solve the dual of a problem and then to deduce the required answer from that solution. 

3.7 Calculation of resistance by energy nnethods 

It has already been observed that capacitances can be calculated by energy methods to an accuracy 
which is adequate for many purposes. Similar methods can be used to calculate resistances. In this 
case it is assumed that the effect of perturbing the distribution of the electric field and the current 
density from equilibrium is to increase the rate of energy dissipation. Now the power dissipated in a 
resistance can be written as 

W = PR = V^/R (Equation (3.7)) 

so, if the total current is given, and the power dissipated is calculated as W from an approximate 
distribution of the current, then 

W'>I'R so that R<W'/l^ (3.19) 

giving an upper bound for the resistance. If, on the other hand, the potential difference across the 
resistance is given and the power is calculated from an approximate set of equipotentials then 

W">V^Ir so that R>V^IW" (3.20) 

giving a lower bound for R. 

3.8 Sunnnnary 

In this chapter we have explored the links between the theory of electromagnetism and the theory of 
electric circuits. It has been shown that Ohm's law and Kirchhoff s laws are special cases which arise 
when electromagnetic theory is applied to lumped components and to electrical networks. 

On occasion it is necessary to calculate the resistance of a piece of conducting material. Three 
possible approaches have been discussed: direct calculation, the use of Laplace's equation and the use 
of energy methods. Direct calculation is possible only in a limited number of cases where the current 
density is uniform over surfaces of simple shape. The application of Laplace's equation can, in 
principle, solve all problems involving uniform conducting materials. Similar numerical methods can 
be used when the materials are not uniform. These methods are available as computer packages of 
great sophistication. For many purposes great accuracy is not required and then the energy methods 
are of value because of their ability to provide acceptable accuracy with little effort. 
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4. The magnetic effects of electric currents 

Objectives 

• To explain the relationship between magnetism and electricity. 

• To introduce the Biot-Savart law and the magnetic circuit law and to show how they are used 
to calculate the magnetic flux density produced by electric currents. 

• To introduce the concept of magnetic scalar potential and to discuss the types of problem in 
which it can be used. 

• To show that forces are exerted on current-carrying conductors in magnetic fields. 

4.1 The law of force between two moving charges 

The theory of magnetism was developed independently from that of electricity until Oersted (1820) 
showed that the two subjects were linked. This link was later explored by Ampere, Faraday and 
Maxwell, but it was not until the theory of relativity was developed by Einstein that the relationship 
between them could be fully understood. It is now clear that magnetic effects can be regarded as a 
consequence of the motion of electric charges. It is not necessary to postulate the existence of 
magnetic poles or dipoles to explain experimental observations. The magnetic properties of materials, 
discussed in Chapter 5, can be explained by assuming the existence of circulating currents on the 
atomic scale. In this chapter we shall confine our attention to magnetic effects in free space. 
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The figure on the right shows two charges which have 
constant velocities vi and V2. It can be shown that the force 
exerted on Q2 by Qu is given by 




F = -^r + ^^(v,A(v,Ar)) = F,+F„ (4.1) 

where r is a unit vector pointing from 1 to 2 and the symbol a represents the vector cross product. 
The first term in this equation is the electrostatic force (Equation (1.1)). The second term, which we 
shall recognize as the magnetic force, is a result of the motion of the particles. The constant jU^ is 

known as the primary magnetic constant; it has the value 4n x 10'^ H m \ (The Henry (H) is the unit 
of inductance - see Chapter 6). 

It is instructive to compare the magnitudes of the two forces Fe and Fm. To simplify matters let us 
suppose that Vi is parallel to V2. This gets rid of the vector products without affecting the order of 
magnitude of the result. The ratio of the magnitude of the magnetic force to that of the electrostatic 
force is then 



c 



^o/^oViV2 = ^ (4.2) 



It can be shown that the constant c is the velocity of light (see Chapter 7), so c' has a numerical value 
which is approximately 10^^ m^ s"^. Thus the magnetic force is negligible compared with the 
electrostatic force unless either: 

1 . both particles have velocities close to that of light, or 

2. the electrostatic component of the force is cancelled by the electrostatic force of an equal and 
opposite stationary charge. 

The latter is just what happens when electrons are moving along conducting wires because those 
contain equal numbers of fixed positive charges. Although the electrons move at velocities much less 
than that of light, there are enormous numbers of them, so useful forces are produced. That is why 
magnetic forces can be used to generate electric power and to drive electric motors. Note that the 
cross product in Equation (4.1) means that the magnetic force on Q2 is always at right angles to V2. 

4.2 Magnetic flux density 

In the development of the theory of electrostatics in Chapter 1 the sources of electrostatic force were 
separated from their effects by introducing the concept of the electric field. Following the same 
approach, we introduce a new vector B and split the magnetic term F^ in Equation (4.1) into two 
parts: 
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B = 74KAr) (4.3) 



and 

F™=e2KAB) (4.4) 

The vector B is referred to, somewhat confusingly, as the magnetic flux density. The units of B are 
tesla (T) or webers per square metre (Wb m"^). Equation (4.3) provides the basis for calculating the 
magnetic flux density produced by any combination of charges moving in free space. The effect of 
that flux density on any other charges can then be found by using Equation (4.4). 

For the moment we will concentrate on ways of computing magnetic flux densities. This information 
might be needed to calculate: 

the forces on moving charges; 
the forces on current-carrying conductors; 
the forces on iron surfaces; 
self- and mutual inductances; 
induced voltages and currents; 
eddy current losses. 

The second and third of these are mainly applied to electric machines, and we shall refer to them only 
briefly. The last two are dealt with in chapter 6. 

The source of the magnetic flux is often a current flowing in 
a wire, so another form of Equation (4.3) is needed. 
Consider an element of length d/ in a wire carrying a current 
/ as shown in the figure on the right. Let the current consist 
of charges of magnitude q moving with mean drift velocity 
vj. If there are n charge carriers per unit volume and the 
cross-sectional area of the wire is A, then the current in the 
wire is given by / = nqAv^ from Equation (3.2). But the 

charge contained in the element of wire is nqA dl , so that 

Qv^ = {nqA dl) =Idl (4.5) 

The contribution to the magnetic flux density from the charges within the element is therefore, from 
Equation (4.3), 




dB 



= ^^^^^ (4.6) 



where dl is a vector element of length. 
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But isolated current elements cannot exist so Equation (4.6) must be integrated around the complete 
circuit in which the current is flowing. The result is 



This equation is known as the Biot-Savart law. It can be used to compute the magnetic flux density at 
any point produced by known currents flowing in a given arrangement of conductors. Since the flux 
density is proportional to the current such problems are linear and the principle of superposition can 
be used. 

Figure 4.1 shows a short straight section of wire carrying a current/. To find an expression for the 
magnetic flux density at the point P at a perpendicular distance R from the wire we consider a small 
element of the wire of length dl. 




(4.7) 
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Fig. 4.1 Calculation of the magnetic flux density produced at the point P by a short length of wire 
carrying current /. 



From Equation (4.6) the magnetic flux density at P produced by the element dl is 
dB = 



jUqIcII cos 0 



4n:R' sec' 0 



(4.8) 



But 



dl = d{Rtan0)=Rsec^ Odd 
so, substituting for dl in Equation (4.8), 
jUqI COS 0d0 



dB = 



4nR 



The magnitude of the flux density at P produced by the whole wire is found by integrating over the 
length of the wire 



B = 



AnR J " 



cos0d0 = 



jUqI sin a 
InR 



(4.9) 



Examination of the vector cross product in Equation (4.6) shows that B is perpendicular to both R and 
to the wire, and that it is directed into the paper. 
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4.3 The magnetic circuit law 

A useful method for calculating the magnetic flux density can be derived by considering the line 
integral of the vector B around a closed path. Figure 4.2 shows such a path encircling a long straight 
wire carrying current /. The return wire for the current is assumed to be far enough away for it to have 
no appreciable field anywhere on the path. The flux density at any point on the path can be calculated 
from the Biot-Savart law by allowing a to tend to nil in Equation (4.9), with the result 




Fig. 4.2 A general path of integration around a long straight wire. 

From symmetry, the vector B must be perpendicular to the wire and to the line OP. For a small 
movement dl along the path the contribution to the line integral is 



B'dl = BRdd = ^dd 
In 



(4.10) 



This expression is independent of the radius and so, for any closed path encircling the wire, 



B dl = //o/ 



(4.11) 



This result is known as the magnetic circuit law. It may be stated in words as follows: The line 
integral of the vector B around a closed path is equal to juq times the total current enclosed. This form 
of the law is derived from Equation (4.1 1) by using the principle of superposition. An equation in the 
form of (4.11) is written for the magnetic flux density due to each conductor and then the set of 
equations is summed. The result remains true when the current is distributed as a current density 
rather than being confined in a set of wires. 

s 

The general form of the law can be written 



B dl = //Jj J dA 




(4.12) 
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where the normal component of the current density J is integrated over an open surface S which spans 
the closed curve as illustrated by the figure on the right. This equation, like Gauss' theorem in 
electrostatics, can be applied usefully to problems only when it is possible to deduce the distribution 
of the magnetic flux density from the symmetry of the problem. We shall see in Chapter 5 that it has a 
very important use in the solution of problems in which the greater part of the magnetic flux is carried 
by an iron circuit. 

4.4 Magnetic scalar potential 

Although the theory of magnetism has been developed here without using the idea of magnetic poles, 
it is useful to introduce the magnetic analogue of electric potential. This concept is a convenient 
mathematical device which does not have any physical significance. The magnetic scalar potential 
difference between the points P and Q, by analogy with Equation (1.13), is 



where the constant juq has been introduced so that Equation (4.13) has the same form as the equations 
given in older texts which assume the existence of magnetic poles. 




(4.13) 
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The magnetic scalar potential is of limited usefulness 
because it can be given a unique value only in a 
region of space in which there are no electric 
currents. This point is illustrated by the figure on the 
right. It can be shown from Equation (4.10) that the 
difference in the scalar potential between P and Q is 
independent of the path of integration chosen 
provided that it does not encircle any electric current. 
If, however, a path such as 2 were chosen for the 
integration, the result would be to subtract / from the 
value of the potential difference calculated. 

What is worse still is that, by making several loops around the current, it is possible to make the 
potential difference take a whole series of different values, because the value of the integral changes 
by / for each complete loop made around the current. Thus the magnetic scalar potential can be given 
a unique value only in a region of space containing no currents. Despite this limitation there are still 
times when the idea is useful. 

The differential form of Equation (4.13) is 

B = -//,grad* =-//oV* (4.14) 

Moreover, since it has been assumed that free magnetic poles do not exist, the magnetic analogue of 
Gauss' theorem is 

(^B-dA = 0 (4.15) 

or, in words, the flux ofB out of a closed surface is zero. The unit of magnetic flux is the weber (Wb). 
The differential form of Equation (4.15) is 

divB = V-B = 0 (4.16) 

Combining Equations (4.14) and (4.16), we find that must obey Laplace's equation in free space 
just as the electric potential does. That is 

div(grad* ) = V'* =0 (4.17) 

The commonest use of the magnetic scalar potential is for calculating the magnetic flux density in the 
space between iron pole pieces. We shall return to this subject in Chapter 5. 
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4.5 Forces on current-carrying conductors 

We have already seen that an element of a current-carrying conductor can be regarded as equivalent to 
a moving charge. Substitution from Equation (4.5) into Equation (4.4) gives the force on a current 
element as 



The total force on a circuit carrying an electric current can be found by integrating Equation (4.18) 
over all the current elements in the circuit. This is of central importance in electric power engineering 
but rarely needed in electronics. You should consult other texts for details of its application. 

4.6 Sunnnnary 

The law of force between two electric charges moving with constant velocities differs from the 
electrostatic law of force by the presence of an additional term. This additional component of force is 
the magnetic force. Magnetism is thus shown to be a particular aspect of electricity rather than a 
separate phenomenon. The concept of magnetic flux density is a useful tool for calculating the 
magnetic forces on charges and currents. There are two ways of calculating the magnetic flux density 
produced by electric currents: the Biot-Savart law; and the magnetic circuit law. 

The concept of scalar potential, which is useful in electrostatics, is found to be less useful in 
magnetostatics. It can be given a unique definition only in a region of space which is free from 
current-carrying conductors. The main application of the magnetic scalar potential is in calculating the 
magnetic flux density distribution in air gaps in magnetic circuits. The question of the magnetic force 
on a current-carrying conductor has been touched on only briefly because of its relative lack of 
importance in electronics. 




(4.18) 
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5. The magnetic effects of iron 

Objectives 

• To show that the polarization of iron in a magnetic field can be represented by an equivalent 
surface current. 

• To introduce the magnetic field vector H as a means of calculating magnetic fields in the 
presence of iron without needing to calculate the equivalent surface current distribution. 

• To introduce a revised form of the magnetic circuit law which is valid in the presence of 
magnetic materials. 

• To derive the boundary conditions at the interface between two materials with different 
magnetic properties. 

• To show that iron can be regarded as a conductor of magnetic flux and to consider the 
application of the idea to magnetic screening. 

• To discuss the calculation of fields produced by current-carrying conductors in the presence 
of iron by the method of images. 

• To apply the analogy between electric current and magnetic flux to magnetic circuits, 
introducing the concepts of magnetomotive force and reluctance. 

• To consider approximate ways of dealing with fringing and leakage fluxes in magnetic 
circuits. 

• To discuss the description of the properties of a magnetic material by its hysteresis loop, 
including the concepts of remanence and coercive force and the distinction between hard and 
soft magnetic materials. 
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• To consider the solution of magnetic circuit problems when non-linear effects are important. 

• To consider the solution of magnetic circuit problems involving permanent magnets, 
including questions of the stabilization of the magnets and Evershed's criterion for optimizing 
the design of a circuit. 

5.1 Introduction 

All engineering students must, at some time, have played with magnets and been fascinated by their 
ability to pick up, not just single iron objects, but whole strings of pins or paperclips. The strange 
property of iron and one or two other materials, which we call ferromagnetism, is employed in many 
ways in electronics. 

• The ability of ferromagnetic materials to concentrate and conduct magnetic flux is used in 
inductors and recording heads. 

• Permanent magnets find application in loudspeakers. 

• Special high-resistivity materials known as ferrites are used to make aerials, inductors for 
radios and microwave components such as circulators and isolators. 

• Magnetic media are used for recording purposes in the form of audio and video tapes and 
computer disks though this is less important than it used to be. 

• The most important uses of ferromagnetic materials, if sheer weight of material were to be the 
criterion, are in electrical power engineering, in transformers, generators and motors. These 
lie outside the scope of this book, but they are important examples of the application of the 
fundamental ideas in electromagnetism discussed in this chapter. 

• 

5.2 Ferronnagnetic nnaterials 

The overwhelming majority of materials which exist in the world have magnetic properties which 
differ so little from those of a vacuum that, for nearly all engineering purposes, they are 
indistinguishable from it. The exceptions to this rule are the elements iron, cobalt and nickel, together 
with certain materials which include them. Their special properties are illustrated by Fig. 5.1. 

Figure 5.1(a) shows an air-cored solenoid carrying a current/. The flux density in the vicinity of the 
solenoid, somewhat as shown, could be mapped with the aid of a flux meter. It is worth emphasizing 
that there would be no measurable change in the field if the solenoid were filled with a core made of 
brass, aluminium, wood, plastic, or any other non-ferromagnetic materials. If, on the other hand, an 
iron core were inserted the result would be a considerable increase in the flux density, as shown in 
Fig. 5.1(b). 
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fa) 



(b) 



Fig. 5.1 The nnagnetic fields of (a) an air-cored solenoid and (b) a sinnilar solenoid with an iron core, 
showing diagrannnnatically the great increase in the strength of the field in the second case. 



This property of iron arises because the iron atoms are themselves tiny 
permanent sources of magnetic flux like current loops. If each of these 'loops' is 
represented by an arrow to show its polarity, then a qualitative explanation of 
ferromagnetism can be given. In an un-magnetized piece of iron the atomic 
magnets are arranged, head to tail, in closed loops. All the flux then circulates 
within the iron and none of it outside. The effect of applying an external flux 
density to the iron is to cause the atomic magnets to line up with the external flux 
so that their flux is added to it. These two situations are illustrated 
diagrammatically in the figures on the right. The actual processes of 
magnetization are more complex than this, but there is no room in this book to 
enter into them in more detail. The macroscopic phenomena of ferromagnetism 
are discussed later in this chapter. 



DO 



Ur>magnetized 



Magnetized 



From Fig. 5.1(b) we observe that the field outside the solenoid could just as well be produced by the 
addition of a second winding, with a suitable current and distribution of turns, to the air-cored 
solenoid of Fig. 5.1(a). This suggests that, just as the polarization of dielectric materials can be 
represented by a surface-charge distribution, so the polarization of ferromagnetic materials can be 
represented by a surface-current distribution. This is consistent with the view adopted in this book that 
magnetic phenomena are best regarded as resulting from the motion of electric charges. 
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Mean 
radius fl 



Fig. 5.2 An iron ring of nnean radius R with N turns of wire wound unifornnly upon it. 

To illustrate the effect of the magnetization of iron we consider an iron ring with turns of wire 
wound uniformly upon it as shown in Fig. 5.2. Let the current in the winding be /, and let the 
magnetization of the iron be represented by an additional current /' in the winding. Then, by the 
magnetic circuit law (Equation (4.1 1)) the flux density in the ring is obtained from 



ittrb = ju,n{i + r) 



(5.1) 



or 



B = 



N 



//o (/ + /') 



27rR 



(5.2) 
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Now experiments with electromagnets show that the magnetic flux density produced increases with 
the current, though not necessarily in a linear fashion, so let us assume that f = Xm^ where Xm 
the magnetic susceptibility. 

Then 

B = ^mA^ + X^) (5.3) 



or, if /u^=\ + x„, 

N N 



where the quantities ju and ju^ , are known as the permeability and the relative permeability of the 

iron, respectively. The analogy with the case of dielectric permittivity is only superficial because it is 
hardly ever valid to regard the magnetic permeability as a constant. The other distinguishing feature is 
that the relative permeability of iron is much greater than the relative permittivity of common 
materials, being of the order of lO"^. 

In problems involving dielectric materials we found it useful to make use of the vector D to avoid 
having to calculate the distribution of polarization charges. The same approach can be used in 
magnetic problems. We define a new vector H by the equation 

H = - (5 5) 

so that, for the simple case shown in Fig. 5.2, we have 

// = ^ ,5.6) 

iTlR 

The vector H is known as the magnetic field. It is measured in A m"^ in SI units, as can be seen from 
Equation (5.6). Since the strength of the field is proportional to the product of the current in the 
winding and the number of turns, it is common to speak of a field strength of so many 'ampere-turns 
per metre'. Equation (5.6) also shows that the vector H depends only upon the distribution of electric 
currents and not upon the arrangement of any iron present. 

Making use of the vector H, we can deduce a form of the magnetic circuit law (Equation (4.12)) 
which is valid when magnetic materials are present: 
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|Hdl = JJjdA 

To find the distribution of magnetic flux density in problems involving 
iron we first calculate H from the known distribution of currents, using 
Equation (5.7), and then calculate B in each region from Equation (5.5) 

Magnetic problems cannot usually be solved by the application of potential theory and the use of 
Laplace's equation because of the non-linearity of magnetic materials. For the same reason it is not 
possible to make use of the principle of superposition unless linearizing assumptions can be made. 
This contrasts with the situation in electrostatics where the distribution of charges is frequently 
unknown but the potentials are given on the conducting boundaries. 

5.3 Boundary conditions 

Whenever a magnetic problem involves two or more regions of space with different magnetic 
properties it is necessary to use the relationships between the magnetic vectors on the two sides of 
each boundary. 

The figure on the right shows the boundary between two materials having permeabilities //^ and . 

Consider a Gaussian surface in the form of a 'pillbox' as shown. If the cross-sectional area of the box 
in the plane of the surface is dA, then the flux of B out of the box is 




{B„,-BjdA = 0 (5.8) 

from Equation (4.15), since the contribution of the flux through the sides of the box becomes 
negligible as the height of the box tends to zero. The consequence of Equation (5.8) may be stated in 
words: The normal component ofB is continuous across a boundary. 

The boundary condition on H is found by applying the magnetic circuit law to 
a small closed path encircling part of the boundary as shown in the figure on 
the right. If the length of the path parallel to the boundary is dl, then the line 
integral of H around it becomes 

{H^^ - )dl = current enclosed (5.9) 

It should be noted that the current referred to is an actual electric current. It does not include the 
hypothetical surface currents which were introduced as a convenient way of describing magnetic 
materials. In many cases the surface current will be zero and then the boundary condition is that the 
tangential component ofHis continuous across the boundary. 




74 



Electromagnetism for Electronic Engineers 



5. The magnetic effects of iron 



5.4 Flux conduction and nnagnetic screening 

Equation (5.9) has one very important practical corollary. Consider an iron bar surrounded by air. 
Assuming that there is no surface current flowing along the bar, we know that the tangential 
components of H are equal inside and outside the bar. But, since B = //H , the component of the 
magnetic flux density parallel to the bar is much greater inside it than outside it. 

If such a bar is placed in a magnetic field in which the total flux is 
fixed, with the axis of the bar aligned with the field, then the 
requirement that Biron » Bair means that most of the fiux will pass 
through the iron as shown in the figure on the right. The flux density 
in the air is then much less than it was before the iron was introduced. 
The iron therefore acts as quite a good conductor of magnetic flux. 

The analogy between flux conduction by iron and electric conduction is strengthened still further if we 
consider what happens to the direction of a flux line entering the iron from the air, as shown in Fig. 5.3. The 
boundary conditions require that Bni = B^a and that B^ = iirBta- Now jir, is typically in the range 10^ - 10^, so the 
flux line must turn through an angle even greater than that indicated in Fig. 5.3. The angle of incidence is 
greatest when the flux line within the iron lies nearly parallel to the surface. A simple calculation shows that if 
ai = 89° and jUr= 1000, then a« = 3.3°. For many purposes it is adequate to make the approximation that iron 
is a perfect magnetic conductor, so that all the flux lines enter it at right angles to the surface. 
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Fig. 5.3 The change in direction of a nnagnetic flux line when passing fronn air into iron. 

The ability of iron to conduct flux is put to work in a number of ways. It allows magnetic circuits to 
be constructed which conduct the flux wherever it is required - an idea which we shall consider in 
detail later on. It also allows flux to be excluded from regions where it is not required. Figure 5.4 
shows what happens when a hollow iron cylinder is placed in a magnetic field. The flux lines prefer to 
pass through the iron and the space within the cylinder is almost entirely free from magnetic flux. 
This property can be used to shield sensitive equipment from magnetic fields, and special alloys such 
as mumetal are used for this purpose. 




Fig. 5.4 An iron cylinder acts as a conductor of magnetic flux and screens the volume inside it from 
the external flux. 



The effectiveness of a screen can be expressed in decibels as 



=201ogio 



at a point in the absence of the screen^ 
B at the same point with the screen 



(5.10) 



A mumetal cylinder 100 mm in diameter and 1 mm thick has a magnetic screening effectiveness of 
60dB provided that the external magnetic field is not strong enough to saturate it. It should be noted 
that it is only low-frequency fields which are screened out by the mechanism of flux conduction. A 
different mechanism operates at radio frequencies. 

The fact that iron acts as a good conductor of flux means that it is often possible to calculate the 
distribution of flux around iron pole pieces by regarding them as equipotential surfaces and solving 
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Laplace's equation (4.17). This method can only give relative values of the flux density because the 
potentials of the pole faces are not known. The absolute values can be found by using the fact that the 
flux crossing the air gap is equal to that in the iron, apart from any leakage flux, and applying the 
magnetic circuit law to the whole flux path. The application of the magnetic circuit law to problems of 
this type is discussed in the next section. 

The fact that an iron surface can be regarded as a magnetic equipotential surface allows the method of 
images to be used for the solution of problems in magnetostatics. 

5.5 Magnetic circuits 

The analogy between the conduction of electricity by conductors and the conduction of magnetic flux 
by iron leads us to the useful concept of a magnetic circuit. Figure 5.5 shows a simple magnetic 
circuit formed by a square iron core with a narrow air gap in it at ^. A coil of turns of wire is wound 
on the core and carries current /. To make the problem easy to handle we make some simplifying 
assumptions. Let us suppose that the magnetic field strength is the same everywhere within the iron, 
having the value Hi, and that the magnetic field in the air gap is Ha. 



Air gap 

Fig. 5.5 A simple magnetic circuit made up of an iron core with an air gap in it. The flux is supplied by 
a winding of N turns of wire. 

Now consider a closed path around the circuit which follows the centre line (shown dotted). Applying 
the magnetic circuit law to this path gives 



Strictly speaking, the length of the iron path is {AL - g) but, provided that g is small, the difference is 
negligible given the other assumptions which have been made. 

The boundary condition (5.8) requires that 




4LH,+gH^=NI 



(5.11) 




(5.12) 
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where O is the total flux circulating, and A is the cross-sectional area of the iron at right angles to the 
direction of the magnetic field. We have therefore assumed that all the flux due to the coil is contained 
within the bounds of the iron core and its projection across the air gap. We also know from Equation 
(5.5) that 



B^=MH, and B^=MoH^ 



(5.13) 



Substitution for //, and Ha in Equation (5. 1 1) in terms of O and the various constants gives 



(5.14) 



which can be written 



where -|| (= NI ) is known as the magneto-motive force and as the reluctance. 



(5.15) 
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Equation (5.15) is analogous to Ohm's law for electric circuits, but the analogy must not be pressed 
too far. Unlike an electric current, the magnetic flux in the circuit dissipates no energy and there are 
no circulating magnetic charges. Moreover, it cannot be emphasized too often that magnetic materials 
do not normally behave in a linear fashion, so the reluctance of a circuit can only be regarded as 
approximately constant at best, and that only under a limited range of conditions. We will return to 
give this point further consideration after discussing the behaviour of real magnetic materials. First, 
though, it is useful to mention two other matters: fringing and leakage of magnetic flux. 



5.6 Fringing and leakage 



In the discussion of magnetic circuits above it was assumed that the flux lines 
in the air gap passed straight across, as shown in the upper figure on the right. 
A little thought shows that this is an over simplification and that the fiux 
lines must actually spread out into a fringing field, as shown in the lower 
figure in the margin. The truth of this observation can be demonstrated by 
laying a sheet of paper over the air gap and sprinkling iron filings on it. 
Clearly the reluctance of the magnetic circuit is affected by the fringing of 
the field around the gap so we need an estimate of how big the effect is. 











//// 




N 



A detailed solution of the problem would require the solution of Laplace's equation in the air gap 
(assuming the iron surfaces to be equipotentials). This is unnecessarily complicated for most 
purposes. A usefiil rule of thumb is that the effective cross-sectional area of the air gap is larger than 
that of the pole face by a strip of width g/2 around it, as shown in Fig. 5.6. Equation (5.12) then 
becomes 



and the reluctance of the circuit is 



4^. + . ^ 



(5.16) 



(5 17) 




Fig. 5.6 The fringing of the magnetic field in an air gap can be allowed for, approximately, by adding 
an extra strip around the area of the pole face. 
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Fringing fields are put to use in magnetic recording heads. Figure 5.7 shows a typical arrangement. 
Most of the flux in the magnetic circuit passes directly across the air gap, but some of the fringing 
field passes through the magnetic coating on the tape and magnetizes it. To achieve linear 
performance the signal supplied to the coil is superimposed on an a.c. bias current whose frequency is 
at least five times the maximum signal frequency. 

The tape is replayed by passing it across either the same head or a similar read head where the varying 
magnetization of the tape produces a varying flux in the magnetic circuit. We shall see in Chapter 6 
that varying the flux passing through a coil causes an electromotive force to be induced in it. The 
e.m.f. can be processed to recover the original signal waveform. 




Tape Y/////////////////A ^ 

\ 

Tape backing 

Fig. 5.7 The arrangement of a magnetic recording head. 

A second difficulty in calculating the flux in the air gap is illustrated by Fig. 5.8. Some of the flux 
simply takes a short cut such as A- A' and never reaches the gap at all. This flux is known as leakage 
flux. It exists because, although iron conducts flux better than air, the air is a poor magnetic insulator. 




Fig. 5.8 An electromagnet showing a typical flux leakage path A-A\ 

If we wish to pursue the analogy with electric circuits it is necessary to think of the material 
surrounding the conductor as being a poor conductor (such as wet earth) rather than an insulator. 
Fortunately, in most cases, a very crude estimate of the leakage flux is adequate to provide a 
satisfactory value for the flux in the air gap. 
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5.7 Hysteresis 



In the preceding sections the permeability of iron has been taken to be a constant. While this is a 
useful first approximation, which can be justified in a limited number of cases, it is not generally 
possible to regard permeability as a constant. The relationship between the magneto-motive force 
applied to a specimen of iron and the magnetic flux density produced is a far from simple one. The 
magnetization of the iron depends not just on the present value of the magneto-motive force, but also 
on the previous history of the specimen. In this book we are restricting our attention to macroscopic 
phenomena, so no attempt is made to offer an explanation of their causes in terms of the atomic theory 
of matter. 

One way of investigating the relationship between the magneto-motive force applied to a specimen of 
iron and the flux density produced would be to use a magnetic circuit like the one shown in Fig. 5.5. 
Provided that the air gap is very small, it is possible to neglect the effects of fringing. The flux density 
in the air gap is then the same as that in the iron and its strength could be measured by placing a Hall 
probe in the gap. The total magneto-motive force applied to the circuit is M, so the magneto-motive 
force applied to the iron is 




(5.18) 
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Now -|| / is the product of Hi and the length of the part of the circuit lying within the iron, so it 
would be possible to deduce values of Hi from measurements of the current in the coil and the flux 
density in the air gap. The results of such measurements are customarily displayed by plotting Bi 
against Hi. There are better ways of obtaining the information than the experiment described, but we 
are concerned here with the properties of the iron rather than the details of how they are measured. 




Fig. 5.9 A typical hysteresis loop showing the initial magnetisation path (at)), saturation (at b), 
remanence (at c), the coercive force (at d), and an inner loop (shown dotted). 

Figure 5.9 shows a typical 5-// curve. If the specimen of iron is initially un-magnetized, then its state 
can be represented by the point a at the origin of the graph. Now suppose that the current in the coil is 
gradually increased. It is found that the B-H curve follows the path ab, rising slowly at first, then more 
rapidly, and, finally, levelling off This curve is known as the initial magnetization curve of the iron. 
When the curve levels off at high values of Hi the material is said to be saturated. 

If the current in the coil is reduced to zero, the flux density does not retrace the curve ab. Instead, it 
follows a curve such as be, so that the flux density in the circuit does not fall to zero when the current 
is zero. The value of Bi at c is referred to as the remanence of the iron. 

In order to reduce the flux in the circuit to zero it is necessary to reverse the current in the coil. The 
value of Hi needed to make Bi zero (the point d) is known as the coercive force. Increasing the current 
in the coil beyond the level needed to reach d eventually produces saturation of the iron at e with the 
direction of the flux opposite to that at b. Finally, reducing the current to zero and then increasing it 
with reversed polarity produces the curve efgb. 

If the current in the coil is repeatedly taken through the same cycle, the hysteresis loop bcdefgb is 
traversed repeatedly in a stable manner, provided only that the maximum value of the current in the 
coil is the same for both polarities. In many applications the flux in the iron is produced by a periodic 
(though not necessarily sinusoidal) current, and the iron then behaves in the manner described. If the 
maximum current is increased still further, there is no change to the loop because the magnetization of 
the iron cannot be increased beyond its value at saturation. If, on the other hand, a periodic current of 
smaller amplitude is used, it is found that the behaviour of the iron is described by a smaller loop such 
as the one shown dotted in Fig. 5.9. 
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The properties of magnetic materials can be varied by making alloys with different proportions. They 
are also strongly influenced by the ways in which the finished material is prepared, especially any 
heat treatments used. A very great variety of materials now exists. Their properties can be looked up 
in reference books. From a practical point of view the main division among magnetic materials is 
between 'hard' and 'soft' materials. The differences between these two classes are illustrated in Fig. 
5.10. 




H 



Fig. 5.10 Typical hysteresis loops for hard and soft magnetic materials. 

Soft magnetic materials are characterized by having narrow hysteresis loops, low remanence and 
small coercive forces. They are therefore easily magnetized and demagnetized and are used as 
conductors of flux in magnetic circuits and magnetic screens. Because these materials have narrow 
hysteresis loops their behaviour can be approximated by their initial magnetization curves. Moreover, 
if the loop is also nearly straight for flelds below saturation, like the one shown in Fig. 5.9, it is 
possible to make the approximation that the permeability is constant. 

Hard magnetic materials have broad hysteresis loops, high remanence and high coercive forces. 
Such materials are difflcult to demagnetize. They are used for making permanent magnets and 
magnetic recording materials. 

A further important feature of the properties of magnetic materials is illustrated by Fig. 5.1 1. Suppose 
that the material has been magnetized to saturation {b) and then demagnetized to the point p by the 
application of a reverse fleld. If the magnetizing fleld is then reduced to zero, the working point of the 
material moves to q along a minor loop. Minor loops are usually narrow and it is possible to 
approximate the loop between p and q by the straight line pq, which is usually nearly parallel to be. 
Provided that the demagnetizing fleld applied does not take the working point past p back onto the 
main loop, the material will operate in a stable manner along the minor loop pq. This fact is of 
considerable importance in the application of permanent magnets. 
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Fig. 5.1 1 Part of a B-H curve showing a nninor loop p-q. 

5.8 Solution of problems in which /j cannot be regarded as constant 

Consider again the magnetic circuit shown in Fig. 5.5. Neglecting leakage and fringing fluxes we may 
assume that 5, = 5^. The application of the magnetic circuit law to the circuit gives 

NI = 4LH. + gH^ = 4LH, + & 

Mo 

or 

B,=^{NI-4LH.) (5.19) 
g 
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A second condition is that Bi and ///must be related to each other by the hysteresis loop of the 
material: 

B,=f{H,) (5.20) 

Equations (5.19) and (5.20) are a pair of simultaneous equations in 5/ and Hi. However, because (5.20) 
is not the equation of a straight line, they are non-linear simultaneous equations. 

They can be solved graphically, or by using an analytical approximation to the equation of the 
hysteresis loop. The graphical solution is shown in Fig. 5.12. The working point Q is at the 
intersection of the straight line represented by Equation (5.19) with the hysteresis loop. The intercept 
of this line on the horizontal axis is the total magneto-motive force in the circuit, while its slope is 
determined by the relative sizes of the iron path and the air gap. The method is similar to the process 
of finding the working point of a transistor by drawing a load line across the characteristic curves. 




4L 

Fig. 5.12 A magnetic circuit problem can be solved by drawing a load line across the B-H plot for the 
iron. The working point of the circuit is at Q. 

5.9 Permanent magnets 

We have already noted that some magnetic materials can be magnetized so that they produce a 
substantial magnetic flux even when the magnetizing field is removed. These hard magnetic materials 
are used for making permanent magnets. Permanent magnets are used in loudspeakers and in a 
variety of microwave devices such as ferrite isolators, magnetron oscillators and travelling-wave tube 
amplifiers. 

Figure 5.13 shows a typical B-H plot for a permanent magnet material. This curve, which is the part 
of the hysteresis loop lying in the second quadrant, is known as the demagnetization curve of the 
material. 
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H 



Fig. 5.13 The working point for a pernnanent nnagnet circuit lies in the second quadrant of the B-H 
plot at a point such as Q. A pernnanent nnagnet can be stabilized against dennagnetization by external 
fields by operating it at a point Q' on a nninor loop such as R-S. 



Suppose that the circuit shown in Fig. 5.5 is made of this material and magnetized to saturation by 
passing a current through the coil with the air gap bridged by a piece of soft iron. When the current in 
the coil is reduced to zero the working point will be at P in Fig. 5.13. Neglecting fringing and leakage 
fluxes we have, when the soft iron is removed from the gap, 



f— 1 

Vs J 



(5.21) 



from Equation (5.19), setting /= 0. This load line passes through the origin and has negative slope. It 
is plotted in Fig. 5.13 as the line OQ. The effect of opening the gap is therefore to demagnetize the 
magnet to some extent. For this reason it is common to speak of the demagnetizing field of the air 
gap. 

If the piece of soft iron is reinserted in the air gap, the operating point of the magnet moves along a 
minor loop to the point P'. In this state the magnet is largely immune to the effects of external fields. 
For this reason it is usual to store permanent magnets with their air gaps bridged by a piece of soft 
iron known as a 'keeper'. 

Suppose that someone has forgotten to put a keeper on the magnet and that it is exposed to a magnetic 
field which tends to demagnetize it. The working point might then move to R. When the field is 
removed the working point lies at on the minor loop RS. This illustration shows that it is necessary 
to treat permanent magnets carefully if their properties are not to be affected. For example, it is 
unwise to try to force a pair of magnets together with their fields opposing because they are then 
trying to demagnetize each other. 
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It is not a good thing that the working point of the magnet should be sensitive to external influences in 
this way. For this reason permanent magnets are often stabilized by deliberately demagnetizing them 
beyond their working points on the main hysteresis loop. The operation of magnets which have been 
stabilized is stable along a minor loop (or recoil line) such as RS unless the external demagnetizing 
field is very strong. Smaller external fields produce a temporary shift in the working point, but it 
returns to Q' when the perturbing field is removed. 

5.10 Using permanent magnets efficiently 

In a circuit involving a permanent magnet, the cost of the magnet is usually a considerable part of the 
cost of the whole circuit. It is, therefore, desirable to use the magnet material as efficiently as possible. 
A rule of thumb for this purpose can be derived by considering the circuit shown in Fig. 5.5. Let the 
lengths of the magnet and the air gap be // and 4 and let their cross-sectional areas be Ai and Aa. The 
magnetic circuit law gives 



0 



(5.22) 



and the conservation of flux gives 



a 



(5.23) 
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The volume of the magnet is 



V =IA^ =-lA, 



I II 



Now laAa is the volume of the air gap which may be taken as given, and Ha and Ba are the field 
strength and flux density in the air gap. So, the minimum magnet volume to provide a given flux 
density in a given air gap is achieved by making the product HiBi as great as possible. This is known 
as Evershed's criterion. The quantity HiBf is known as the energy product of the magnet for reasons 
which will become apparent in a later chapter. 

Figure 5.14 shows the curves of B against H and BH against B as they are usually plotted. The 
optimum working point of the material is then that which gives the value of B corresponding to the 
maximum value of BH as shown at Qopt. 
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Fig. 5.14 Typical curves of H and BH versus B for a permanent magnet material showing the choice 
of working point at the maximum value of BH to make the most efficient use of the magnet. 



5.11 Sumnnary 



In this chapter we have considered the properties of ferromagnetic materials. We have seen that their 
principal characteristic is that they acquire a very strong magnetization when placed in a magnetic field. As 
with dielectric materials it is possible to represent the properties of the materials by introducing a new vector, 
in this case H, which includes the effects of the polarization of the materials. Using the vector H a new form 
of the magnetic circuit law was derived, which is valid even when magnetic materials are present. 

Consideration of the boundary conditions at the surface of a ferromagnetic material showed that such 
materials can usually be regarded as good conductors of magnetic flux. This property allows the path 
of the flux to be controlled so that it is directed to the place where it is needed, as in a recording head, 
or directed away from sensitive areas in magnetic screening. 

In a limited range of cases it is possible to regard ferromagnetic materials as being linear, but, in 
general, it is necessary to take account of their non-linear behaviour as described by their hysteresis 
loops. 
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6. Electromagnetic induction 

Objectives 

• To show that an electromotive force is induced in a conductor moving through a magnetic 
field, and that a current flows if the circuit is completed. 

• To show that an electromotive force is induced in a loop of wire moving through a non- 
uniform magnetic field and that, by implication, an e.m.f is induced in a loop of wire when 
the magnetic flux linked to it changes. 

• To generalize the first two objectives in the form of Faraday's law of electromagnetic 
induction. 

• To demonstrate the links between the circuit concepts of self and mutual inductance and 
electromagnetic field theory. 

• To discuss the causes of electromagnetic interference and ways of reducing It. 

• To introduce methods, including energy methods, for calculating self and mutual inductances. 

• To introduce the analogy between Z, C and R in two-dimensional problems. 

• To discuss the concept of the storage of energy in a magnetic field, and to show the 
equivalence between calculations of the energy stored in an inductor from field and circuit 
points of view. 

• To discuss the special case of energy storage in a magnetic field in iron and to introduce the 
idea of hysteresis loss. 

• To introduce the induction of eddy currents in a conductor by a changing magnetic field. 

6.1 Introduction 

In this chapter we turn our attention to phenomena involving conductors and magnetic fields in 
relative motion, and magnetic fields which are changing with time. These effects, first investigated by 
Faraday in 1831, complete the link between magnetism and electricity. We have already seen how 
magnetic fields are produced by electric currents. We shall now show that electromotive forces and 
electric currents can be produced by electromagnetic induction. These ideas are fundamental to the 
generation of electricity by electromagnetic machines and to the use of transformers in the distribution 
of electric power. As far as electronic engineers are concerned their importance lies in the fact that 
they provide an explanation of self and mutual inductance, and the means to calculate these circuit 
parameters. 

6.2 The current induced in a conductor nnoving tlirougli a steady 
nnagnetic field 

Consider a piece of straight wire moving with constant velocity v through a 
steady magnetic field B which is at right angles to the wire and to its direction 
of motion, as shown in the figure on the right. The conduction electrons in the 
wire are carried through the field with the wire so, from Equation (4.4), each 
electron experiences a magnetic force qvB acting upon it. The cross product 
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symbols have been omitted because the directions of the wire, the velocity of the wire and the field 
have been chosen to be mutually perpendicular. As a result of this force the electrons are displaced so 
that the ends of the wire acquire polarization charges, as shown in the figure. 

The polarization of the wire produces an electric field within it which opposes the motion of the 
electrons. In equilibrium the electric force on each charge must be equal and opposite to the magnetic 
force, so that 



Now suppose that the ends of the wire are connected to parallel rails by sliding contacts and that the 
circuit is completed by a load resistor, as shown in Fig. 6.1. The potential difference between the ends 
of the rod acts as an electromotive force in the circuit and a steady current flows such that 



E = vB 



(6.1) 



The magnitude of the potential difference between the ends of the wire is therefore 



V = \Edl = vBl 



(6.2) 



IR = vBl 



(6.3) 
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The electric power dissipated in the resistor is 



rR = vBlI 



(6.4) 



The mechanical power input is the product of the force needed to move the wire and its velocity. 
Since the current in the circuit is /, the force needed to move the wire is BlI, from equation (4.18), so 
the mechanical power input is 



Fv^vBlI 



(6.5) 



This is exactly equal to the electric power, demonstrating that this device converts mechanical energy 
into electrical energy. 



C 



B 

® 



J 



Fig. 6.1 A current is induced in a closed circuit in a magnetic field if the shape of the circuit is 
changed. 

Examination of Fig. 6.1 shows that the force exerted on the current-carrying slider by the field is in 
the opposite direction to v. This internal force must be opposed by an equal and opposite external 
force if the wire is in steady motion. The external force is in the same direction as v, so it does work 
on the slider, providing a power input to the device. This simple arrangement forms the starting point 
for the discussion of the generation of electricity by electromagnetic machines. Important as the topic 
is, it is outside the scope of this book and will not be pursued any further here. 

Recalling the definition of magnetic flux which was used in earlier chapters, we can see that, with the 
direction of velocity shown, the flux linked to the circuit of Fig. 6.1 is changing at a rate given by 



JO 
dt 



= -vlB 



(6.6) 



but, from Equation (6.2), this is just the electromotive force developed in the circuit. The direction of 
the electromotive force is such that when a current flows in the circuit the magnetic forces oppose the 
external mechanical force which is causing the flux to change. These ideas can be combined by 
writing 



dt 



(6.7) 
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6.3 The current induced in a loop of wire nnoving through a non- 
unifornn nnagnetic field 

In the experiment described in the previous section it was shown that a current could be induced in a 
circuit by altering the shape of the circuit and so changing the flux linked to the circuit. It is 
interesting to see whether a current could also be induced by keeping the shape of the circuit fixed and 
varying the strength of the magnetic field. This could be achieved by taking a square loop of wire, as 
shown in Fig. 6.2(a), and moving it with a steady velocity v through a magnetic field perpendicular to 
the plane of the loop whose variation with position in the direction of motion is given by the graph of 
Fig. 6.2(b). 





{b) 



Fig. 6.2 A current is induced in a closed loop of wire when it is moved through a non-uniform 
magnetic field at right angles to the plane of the loop. 

At the instant when the wire AD is at x relative to the origin, the e.m.f generated in it is 



^ = vIB 

AD y/^^o 



(6.8) 



where the order of the subscripts shows that the polarity of the e.m.f is such as to drive a current from 
Ato D through the rest of the circuit. In the same way the e.m.f generated in BC is 



^ =vlB 



fx+r\ 
[ I J 



(6.9) 



The remaining sides of the loop are moving parallel to their own directions so no potential differences 
are generated between their ends. The net e.m.f in the circuit is therefore 



■^=vlBr 



x + l X 



V 



/ / 



(6.10) 



directed so as to tend to drive a current in a counter-clockwise direction. Now the flux linked to the 
circuit is 
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f 



/ J 



0 



(6.11) 



V 



so that 



JO 



dx 



= vlB, 



(6.12) 



dt 



0 



dt 



0 



So, in this case also, the e.m.f. is equal to the rate of change of flux linkage. The direction of the flux 
produced by the induced current is opposed to the original field, so that it tends to prevent the flux 
linked to the loop from increasing as the loop is moved. Since it is the relative motion of the loop and 
the field which produces the e.m.f it is reasonable to suppose that the same result would have been 
produced if the loop had been held fixed and the source of the field moved with velocity v in the 
negative x-direction. 

6.4 Faraday's law of electromagnetic induction 

In the previous sections it has been shown that, at least in the two cases considered, an electromotive 
force is induced in a circuit when the magnetic flux linked to the circuit is changed. This phenomenon 
was observed experimentally by Faraday, who generalized his flndings in the law which bears his 
name. 
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Faraday's law of electromagnetic induction states that, if the flux linked to a circuit is changed in any 
way, then an electromotive force is induced in the circuit whose magnitude is proportional to the rate 
of change of the flux linkage to the circuit. 

The term flux linkage is used here to include circuits such as coils in which the wire encircles the 
path of the circuit more than once. In these cases the flux linked to the circuit is the sum of the flux 
linkage to the individual turns. It is often possible to make the approximation that the flux linkage to 
each turn of a coil is the same, so that the total flux linkage is the product of the flux linked to one 
turn and the number of turns. Denoting flux linkage hy A, this relationship is expressed by the 
equation 

A = A^O (6.13) 

We have also seen that the direction of the e.m.f induced in a circuit by a changing magnetic flux 
linkage is always such that it tries to oppose the change of flux which causes it. This statement is 
known as Lenz's law. 

Faraday's law can be expressed in mathematical terms by writing 

-=-^ (6.14) 
dt 

For a simple circuit such as that shown in Fig. 6.3, the flux linkage is given by 

A = jB-dA (6.15) 

s 

where the integration is carried out over an open surface spanning the circuit and the elementary 
vectors dA are taken to point upwards as shown. 




Fig. 6.3 The direction of the current induced in a loop of wire by a magnetic flux which is increasing 
with time. 
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The e.m.f. in the circuit can be calculated in a similar manner by making use of Equation (3.17). The 
direction of the line integral around the circuit is chosen to be in a right-hand corkscrew sense with 
respect to the direction of the vector dA. This direction is shown in Fig. 6.3 by the direction of the 
elementary vector dl. Thus 

^=1 E-dl (6.16) 



Now, if B is in the direction shown and it is increasing in magnitude with time, Lenz's law requires 
the direction of the induced current to be opposite to the direction of the line integral as shown in Fig. 
6.3. 

Substituting for mdA in Equation (6.14) produces a general mathematical form of Faraday's law: 

(f E.dl = -^jjB.dA (6.17) 

Although this expression may seem abstract and complicated it should be remembered that it is no 
more than a restatement of Faraday's law and Lenz's law in the symbolic language of mathematics. 

When a circuit has more than one turn, the integral on the right-hand side of Equation (6.17) is taken 
over all the turns. In the problems which are of interest in electronics the circuit is normally fixed and 
the magnetic field changing with time. Under these circumstances the differentiation can be brought 
inside the integral to give 

fE-dl— n^-dA (6.18) 



6.5 Inductance 

We are now in a position to consider the third of the trio of passive circuit parameters, resistance, 
capacitance and inductance. Figure 6.4 shows two loops of wire arranged so that when a current is 
passed through loop 1 some of the flux produced passes through loop 2 and vice versa. 




Fig. 6.4 Two loops of wire are inductively linked when some or all of the magnetic flux produced by 
one of them passes through the other. 
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Suppose that the flux Hnked to loop 2 is A21 when the current in loop 1 is /i. If the loops are in a region 
of constant permeability, thenyi2i is proportional to /i. The constant of proportionality is 

4i=4^ (6.19) 

where L21 is the mutual inductance coupling the loops. From Equation (6.19) it is evident that 
inductance has the dimensions weber per ampere, and this unit is called the Henry (H). The ratio of 
the flux linked to loop 1 to the current in the same loop is the self-inductance of that loop, given by 

Ai=4^ (6.20) 

Two similar equations can be derived by assuming that the source current is in loop 2. 
These equations can then be brought together by using the principle of superposition to give 

Ai=Aii+Ai2=AiA+A2A 

A2 = + = ^2iA + ^iih (6.21) 
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As this is a passive linear system the matrix of the coefficients of inductance must be symmetrical; 
that is 

Z21 =L^2=M (6-22) 

where Mis the mutual inductance of the loops. Using this notation the double subscripts can be 
dropped and the self inductances represented by the symbols Li and L2. When Faraday's law is applied 
to equations (6.21), we obtain 

^ dL ^ ^ dL 
dt dt 



dl^ dl^ 

2 

dt dt 



= -M^-L^^ (6.23) 



where the positive directions of the currents and the e.m.f s are the same. These are the familiar circuit 
equations for inductively coupled circuits. It is important to remember that these equations hold only 
when the system is linear, that is, when the fluxes are directly proportional to the currents producing 
them. Since the behaviour of iron can be highly non-linear, as we saw in Chapter 5, it cannot be taken 
for granted that iron-cored inductors behave as linear circuit elements. 

The self and mutual inductances in a circuit are related to each other, as will now be demonstrated. 
Consider again the pair of loops shown in Fig. 6.4, but suppose that loop 1 is a coil of Ni turns and 
coil 2 has N2 turns. When a current /i flows in 1 a flux Oi is generated. For the sake of simplicity it 
will be assumed that the whole of this flux is linked to every turn in coil 1, so that 

A, =7ViOi (6.24) 

In general, not all of the flux will pass through coil 2 so that the flux linked to that coil can be written 

A2 =^2^2^i (6.25) 

where is the proportion of the flux generated by coil 1 which is linked to coil 2. It follows 
that ^2 ^1- 

Assuming linearity, we can write 



K = ^ = k^ (6.26) 
If a current is passed through coil 2 we obtain, by a similar argument 

^ = ^ = ^ (6.27) 

Z2 A2 A^2 
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where ki is the fraction of the flux produced by coil 2 which passes through coil 1 . Multiplying 
Equations (6.26) and (6.27) together gives 



(6.28) 



where the coupling coefficient k = sjkjz^ - 1 • 

6.6 Electromagnetic interference 

In Chapter 2 we saw that there can be unwanted capacitive coupling between two circuits. Similar 
problems can arise through inductive coupling. Figure 6.5 shows a simple example of the coupling of 
two circuits through their mutual inductance M. If an alternating current flows in circuit 1, then it 
produces a magnetic flux, some of which may be linked to circuit 2. 




o 



® 



Fig. 6.5 When two circuits are coupled by stray mutual inductance the signal in one circuit can 
interfere with that in the other. For example, if circuit 1 is the a.c. mains a spurious 50 Hz signal could 
be added to the signal \/2 in circuit 2. 

To simplify matters we will assume that I\ is much greater than h so that, approximately. 



and the e.m.f. induced in circuit 2 is 



, 4- dL jo) M / . X 



dt 



(6.29) 



This voltage is added to the signal voltage in circuit 2. The interference is worst when Vi is large and 
the loop impedance of circuit 1 is small. The commonest example of this is when circuit 1 represents 
the a.c. mains and the interference appears as 50 Hz 'mains hum'. In this case it is not possible to 
reduce the interference by altering Vi or (Rsi + Rli) and effort has to be concentrated on reducing M. 
To do this we have to reduce the proportion of the flux produced by circuit 1 which is linked into 
circuit 2. 
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Three strategies are possible: 

1 . Reduce the area of circuit 2 

If the signal source V2 is a transducer some distance from the amplifier, then this can be 
achieved by using a pair of wires twisted together (a 'twisted pair') to connect them together. 

2. Rotate circuit 2 so that its plane is parallel to the flux produced by circuit I 

Rotating one circuit relative to the other can reduce the flux linkage to zero. If the two circuits 
are to be enclosed in the same case in fixed positions this can be a useful approach. 

3. Put a screen of high permeability magnetic material in a position which screens circuit 2 
from the flux of circuit 1. 

The special alloy known as mumetal {jir ~ 10^) is used for this purpose. This method works 
only at low frequencies. 

A particularly troublesome type of electromagnetic interference is caused by earth loops. Figure 6.6 
shows a typical situation in which this problem occurs. Two electronic instruments, A and 5, are 
connected to the main supply by three-core cables whose earth conductors are connected to the 
instrument frames. The instruments are also connected to each other by a coaxial cable whose sheath 
is connected to the frames of the instruments to screen the signal wire from capacitively coupled 
interference. This arrangement produces a closed earth loop as shown by the cross-hatching in Fig. 
6.6. 
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Fig. 6.6 An earth loop, shown shaded, can be fornned by the earth conductors of the electricity supply 
to two electronic instrunnents and the sheath of the coaxial cable connecting thenn together. 

The earth loop has low resistance and large currents can be induced in it by the flux of the mains or of 
some other source of interference. The induced current (//) flows through the resistance of the sheath 
of the coaxial cable (7?c) as shown in Fig. 6.7. This produces a spurious potential difference between P 
and Q which is added to the signal voltage Va- 



Fig. 6.7 The potential difference across the resistance Rc of the sheath of the coaxial cable in Fig. 6.6 
produced by the circulating earth loop current /, is in series with the signal source \/a. 

Two solutions to this problem are possible: 

1 . Break the earth loop so that no current can flow in it. 

This is not as simple as it seems. Disconnecting the earth lead from the plug at i? or 5* is potentially 
very dangerous because lethal voltages could appear on the instrument frames under fault conditions. 
This solution is only possible if a permanent earth connection between A and B can be ensured, for 
example by bolting them to the same rack. Ideally there should only be one path to earth from any 
point on the equipment. The alternative of disconnecting the cable sheath at P means that the signal 
earth follows the path PSRQ. The mutual inductance between this circuit and the mains is high. 

2. Increase the resistance of the earth loop. 

This can be achieved by putting resistors between the signal earth and the frame of each 
instrument as shown in Fig. 6.8. If the added resistors {Rd) are much larger than Rc then the 
circulating current in the earth loop is reduced because the e.m.f is fixed. The unwanted 
voltage appearing across Rc is reduced by a factor which is approximately RcI2Rd. It is 
important that Rd is small enough so that no part of the system which can be touched reaches 
a lethal voltage when the maximum fault current is flowing. \fRc = 0.\ Q and Rd = 50 Q, 
then the interference voltage is reduced by a factor of 1000. If, also, the mains fuse of each 
instrument is rated at 250 mA, then the voltage at P or g cannot exceed 12.5 V with respect to 
earth. 
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Fig. 6.8 The resistance of the earth loop can be increased by connecting resistors Rd between the 
signal earth and the supply earth of electronic instrunnents. 



6.7 Calculation of inductance 



Before proceeding to a discussion of the ways of calculating inductance it is useful to review the 
situations in which such a calculation might be necessary. They are: 

• To find the self and mutual inductances of circuit components such as inductors and 
transformers. 

• To find the inductance per unit length of two-wire transmission lines. 

• To estimate the stray mutual inductance between parts of a circuit such as parallel tracks on a 
printed circuit board. 

Inductance is calculated either by assuming a distribution of currents and computing the flux linkage, 
or by assuming a flux distribution and then finding the ampere turns needed to produce it. The first 
approach is best for many purposes, but we shall see later in the chapter that the two methods lead to 
upper and lower bounds for the inductance when energy methods are used. The calculation of 
inductance in simple cases may be illustrated by two examples. 

Figure 6.9 shows a simple transformer in which the primary and secondary windings having Ni turns 
and N2 turns respectively are wound upon an iron core whose permeability is ju. The mean path length 
of the magnetic circuit is L and its cross-sectional area is A. 
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Let a current /i flow in the primary winding, then, by the magnetic circuit law the magnetic field in 
the core is 

H = ^ (6.30) 



and the flux in the core is 

0 = BA = juHA = t^^l^ (6.3 1) 



The fluxes linked to the primary and secondary windings are 

K,=N,^ = ^I, and A, = iV,0 = (6.32) 

L L 

so that the self and mutual inductances are 



A=^ = ^ and M = Al=/f^VM (6.33) 
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In the same way by assuming that a current h flows in the secondary winding we obtain 

L^^^^t^ and ^^K^ ^N.N^A 

^ L L 

We have shown, incidentally, that Equation (6.22) is valid for this case and that = L^L^ as would 
be expected from Equation (6.28) for the perfect coupling between the windings which has been 
assumed here. 

In real transformers there is always some leakage of flux, so that the inductances calculated in this 
way are only approximate. To achieve greater accuracy it would be necessary to compute the detailed 
distribution of flux in and around the core. The effects of leakage can be minimized by using 
multilayer windings. 

At frequencies above 50 Hz the capacitance between the windings can be important because it limits 
the band of frequencies over which the transformer will work. The complete transformer can be 
represented by an equivalent circuit and methods exist for calculating the parameters of the circuit for 
any particular design of transformer. The core can be regarded as linear to a first approximation if it is 
made of a soft magnetic material and operated well below saturation. Inductors ('chokes') only have a 
single winding. They are sometimes made with an air gap in the magnetic circuit, so that the 
reluctance of the circuit is dominated by the reluctance of the gap and is, therefore, linear. 

Figure 6.10 shows unit length of a coaxial cable whose inner and outer conductors have radii a and b 
respectively. We wish to find an expression for the series inductance per unit length of the line. 




Fig. 6.10 A coaxial cable 



If the ends of the cable are short-circuited then the two conductors form a closed circuit as shown. If 
the magnetic circuit law is applied to a circular path which is concentric with the cable then net 
current is only enclosed when the path lies in the space between the conductors. Thus there is a 
magnetic field only in the space between the conductors and, from considerations of symmetry, it 
must be directed in the tangential direction as shown. 
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Let the current in the conductors be / then using the magnetic circuit law we find that the magnetic 
field strength on a path between the conductors at radius r is 



H = ^ (6.35) 



The conductors form a closed circuit having only a single turn, so the flux linked to it is found by 
calculating the flux crossing a unit length of a radial plane, as shown in Fig. 6.10. Thus, assuming that 
the material in the space between the conductors is non-magnetic 



A = f 



(6.36) 



The inductance per unit length is the flux linkage per unit current 

(6.37) 



In 



The significance of this result will become apparent in Chapter 7 when the use of coaxial cables as 
transmission lines is discussed. The estimation of the inductance per unit length of transmission lines 
of other shapes can be achieved by energy methods. This topic is treated later in this chapter. 

6.8 Energy storage in the nnagnetic field 

In Chapter 2 we saw that the work done in setting up an electric field could be thought of as being 
stored either in a lumped manner in the capacitance of the system ( ^/iCV^) or distributed throughout 
space with an energy density /^D.E. In just the same way the work done in setting up a magnetic field 
can be regarded as being stored either in the circuit inductance or distributed throughout the magnetic 
field. 

It is well known from elementary circuit theory that the energy stored in an inductor is given by 



Jo 9 



2 (6.38) 
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An alternative expression for W is obtained by substituting A/I for L 



W = -AI 
2 



(6.39) 



To derive an expression for the energy density in a magnetic field 
we consider a current loop and a typical flux tube linked to it, as 
shown in the figure on the right. The stored energy associated 
with the tube, obtained fi-om Equation (6.39), is 




AW = -AOI 
2 



(6.40) 



Now, if the cross-sectional area of the flux tube at a particular point is 5A and the flux density at the 
same point is B, then 



(6.41) 
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Applying the magnetic circuit law to the flux tube gives 

I = lll'dl^Y,Hdl (6.42) 



if the flux tube is considered to be made up of volume elements of length dl as shown in the figure 
above. Substituting in Equation (6.40) for AO and /produces 

A r = ^ {BSA)Y, {Hdl) = ^Ya {BHdAdl) (6.43) 

where the vector dot products have been dropped because the vectors concerned are all parallel to 
each other from the definition of the flux tube. {B dA) is a constant of the flux tube, so it is permissible 
to bring it inside the summation. Finally, noting that the product {dAdl) is just the volume, Sv, of the 
element defined by them, the energy stored in the flux tube becomes 

A^ = ^^5//Jv (6.44) 

so that the energy density in the magnetic field is evidently /^BH. As in the case of the electric field 
considered in Chapter 2, a more rigorous derivation shows that the energy density in the field is given 
in general by 

w = ^B-H (6.45) 
so that the total energy stored may be written 

^^IjjjB.HJv (6.46) 
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6.9 Calculation of inductance by energy nnethods 

We have seen that energy methods provide a useful way of obtaining quite good approximations to 
capacitance and resistance in cases which cannot be solved by analytically. These methods can also be 
used to estimate inductances. The argument runs exactly parallel to the cases of capacitance and 
resistance so it will not be given in detail. Starting from Equation (6.38) we discover that 

L<-- (6.47) 



where W is the estimate of the stored energy produced by assuming that the current distribution is 
given. The corresponding lower bound can be obtained by making use of the definition of self- 
inductance to give another expression for the stored energy 

A' 

W = — (6.48) 
2L 

If an estimate W" of the stored energy is obtained by assuming that the flux distribution is given, then 
A' 

L > (6.49) 

6.10 The LCRZ analogy 

From Equation (2.10) it is easy to show that the shunt capacitance per unit length of a coaxial cable is 
'-Wa) 

A comparison between Equations (6.37) and (6.50) shows a marked resemblance. This suggests that, 
if we can calculate the capacitance for a particular two-dimensional arrangement of electrodes, we 
should be able to deduce an inductance for the same geometry from it. The justification for this is that 
the field patterns in each case must satisfy Laplace's equation (1.27). Since this also applies to current 
flowing in a conductor, we can add resistance to the discussion. We will assume that the medium 
between the electrodes has the same properties everywhere. 

It is easy to show that the resistance per unit length for radial current flow between concentric 
electrodes (Fig. 6.1 1(a)) separated by a uniform conducting material of resistivity p is given by 

R = -^\n{bla) (6.51) 
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Fig. 6.11 (a) Coaxial and (b) parallel strip geonnetries 



Similarly, for the parallel strip geometry shown in Fig. 6.1 1(b) we obtain the following parameters per 
unit length if fringing fields are neglected: 



C = sw/d 
L = jud/w 
R = pd/w 



(6.52) 
(6.53) 
(6.54) 
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We deduce that for all such two-dimensional arrangements of electrodes the parameters per unit 
length are given by 

C = sr (6.55) 
L = ju/r (6.56) 
R = p/r (6.57) 

where the parameter T is determined solely by the geometry of the electrodes. For coaxial geometry 

^-W/a) ''''' 

and for parallel strips 

(6.59) 

r = w/d 

By multiplying Equations (6.55) and (6.56) together we obtain 

LC = SjU (6.60) 

which is a constant whose value depends only on the properties of the material surrounding the 
electrodes. 

Dividing Equation (6.56) by Equation (6.55) and taking the square root we get 



Y (6-61) 



where Z has the dimensions of resistance and depends only upon the properties of the material. 

Equations (6.55) to (6.57) can be used to find any two of the three parameters C, L and R for a particular 
geometry once one of them is known. Note carefully that the magnetic flux lines in each case are at right 
angles to the electric field lines. The significance of this will become apparent in Chapter 7. 

These ideas can be extended further by invoking the principle of duality discussed in section 3.6 in 
which the roles of the flux lines and equipotentials are exchanged. It is easy to see from Fig. 6.1 1(b) 
that the parameters of the dual system are given by 

C = s/r (6.62) 
L'^juT (6.63) 
R' = pr (6.64) 



109 



Electromagnetism for Electronic Engineers 



6. Electromagnetic induction 



6.11 Energy storage in iron 

So far our discussion of energy in magnetic fields has been restricted to those cases where the 
permeabiHty is constant. To investigate storage of energy when the permeabihty is not constant, we 
consider an iron ring of mean circumference / and cross-sectional area^ which has turns of wire 
wound upon it, as shown in Fig. 6.12. 




1^ 



Fig. 6.12 The arrangement of a toroidal inductor 

In time dt let the current in the winding increase from Ho {1+ dl) and the flux in the iron increase 
from 0io{0 + d0). The induced e.m.f in the winding is Nd^/dt so the work done in the time 
interval is 



JO 

dW = NI dt = NIdO 

dt 



(6.65) 



But 



and 



Hl^NI 



so 



dW = AlHdB 



(6.66) 



The change in the stored energy when the system is taken from flux density Bi to B2 is given by 



W = Al[^" H dB 



(6.67) 



But Al is just the volume of the iron, so the change of energy density is 



w 



r^2 

= \ HdB 



(6.68) 
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Once again more rigorous argument shows that Equation (6.68) is correct for all cases if the scalar 
product HdB is replaced by a vector dot product H -dB . The work done in magnetizing an initially 
unmagnetized specimen of iron is illustrated by Fig. 6.13. The integral in Equation (6.68) is 
represented by the shaded area between the initial magnetization curve and the vertical axis. 



B 




H 



Fig. 6.13 The work done in magnetizing a piece of iron is represented by the area between the initial 
magnetization curve and the vertical axis. 
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6.12 Hysteresis loss 

A case of particular interest is the work done in taking a piece of iron once around its hysteresis loop. 
In any iron-cored inductor or transformer the core is taken through this cycle for every cycle of the 
a.c. current in the windings. Figure 6.14 shows a typical hysteresis loop. The change in stored energy 
per unit volume in going from 1 to 2 is represented by the shaded area in Fig. 6.15(a). Let us denote 
this energy by wn. 




H 



Fig. 6.14 The work done in taking a piece of iron around its hysteresis loop is proportional to the area 
of the loop. 




Fig. 6.15 (a) The increase in the stored energy of a piece of iron when it is taken around its hysteresis 
loop from 1 to 2; (b) the energy recovered when it is taken from 2 to 3. 

The shaded area in Fig. 6.15(b) represents the change in stored energy per unit volume when the iron 
is taken from 3 to 2, a change which is mathematically possible but physically impossible. This 
energy is W32. The real physical change is obtained by exchanging the limits of integration in Equation 
(6.68) so that W23 = -W32. The change in energy density resulting from going from 1 to 3 via 2 is 

which is represented by the difference between the areas shown in Figs 6.15(a) and 6.15(b), and by the 
shaded area in Fig. 6.14. Since the hysteresis loop is symmetrical the work done in encircling the loop 
once is equal to the area of the loop multiplied by the volume of the iron. It is evident that this work is 
not zero since all real ferromagnetic materials have hysteresis loops with non-zero areas. 
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This phenomenon means that electrical energy is changed into heat in the iron core at a rate 
proportional to the frequency of the signal applied to the device. The conversion of energy in this way 
is known as hysteresis loss. The loss is minimized by using soft magnetic materials for transformer 
cores, as an examination of Fig. 5.10 will show. 



6.13 Eddy currents 



When any conducting circuit is placed in a changing magnetic field currents are induced in the circuit, 
as we saw earlier. This remains true when the circuit is a solid piece of metal. In that case currents, 
known as eddy currents, circulate within the metal, causing power loss by ohmic heating. They are 
put to use in the industrial eddy-current heating (ECH) process, which is used for brazing metal 
components together. In other cases, particularly in transformer cores, they give rise to unwanted 
losses. 

The full theory of the generation of eddy currents is beyond the scope of this book, but it is possible to 
give an approximate treatment of the special case in which the conductor is in the form of a thin strip 
with the magnetic field parallel to its length. 




Fig. 6.16 The element of volume used in estimating the eddy-current loss in a thin conducting strip. 



Figure 6.16 shows a section of a strip of conducting material whose width (/) is much greater than its 
thickness (2t). The strip is placed in a uniform magnetic field directed along the strip with flux density 
5q cos{cDt) . The eddy currents in the strip flow in closed loops like the one shown. From Faraday's 
law we know that the e.m.f induced in the loop is given approximately by 
dB 

^= -21 X — = llxcoB^ ^m{cot) (6.69) 
dt 

since the area of the loop is approximately llx. 



The resistance of the loop is approximately 

R = ^ (6.70) 

dx 

per unit length of the strip, where p is the resistivity of the material. 



The mean power loss in the loop is 
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dW = - 



2R 



X dx 



P 



so that the mean power loss per unit length of the strip is 



W = 



^l(D'Bl~]ct,,_iea/Bl 



P 



1 1^ Jx = 



3p 



(6.71) 



We can also express this as mean power loss per unit surface area 



Ws = 



tWB^ 
3p 



(6.72) 



In deriving this expression we have implicitly assumed that the eddy currents are not strong enough to 
have a significant effect upon the strength of the magnetic field. 
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Equation (6.72) shows that the eddy-current losses depend very strongly on the thickness of the strip. 
This is the reason why transformer cores are made of thin strips or laminations of steel which are 
insulated from each other by a coating of lacquer. They are normally made of special steel which has 
a high resistivity because, as can be seen from Equation (6.72), this also helps to reduce the losses. 
The mean power loss also increases rapidly with increasing frequency. At radio-frequencies the losses 
in laminated iron cores are unacceptably high and magnetic oxides of iron known as ferrites are used 
instead because they have much higher resistivities. 

6.14 Real electronic components 

By this stage it should be apparent to the reader that the behaviour of real electronic components is 
much more complicated than that of the lumped resistors, capacitors and inductors used in circuit 
theory. Whenever current flows through a component a magnetic field is generated, so resistors and 
capacitors must have some self-inductance. The resistance in the winding of an inductor or 
transformer means that there must be some voltage drop across the device even under d.c. conditions. 
The existence of a potential difference between the different parts of the winding implies that there is 
some capacitance present. Although these parasitic effects can usually be neglected at low 
frequencies, they must be included in the circuit representations of the components at high 
frequencies. 

6.15 Sunnnnary 

In this chapter Faraday's law of electromagnetic induction has been introduced. It has been shown 
that an electromotive force can be induced in a circuit which has magnetic flux linked to it either by 
changing the area of the circuit or by changing the strength of the magnetic field. From this starting 
point the idea of inductance was introduced, so establishing another link between electric circuit 
parameters and electromagnetic field theory. 

Stray inductances between parts of electronic circuits can produce unwanted coupling. The causes of 
electromagnetic interference were considered, and some cures suggested. The calculation of 
inductance was discussed using both direct calculation and energy methods. The latter are based on 
the idea of the stored energy associated with a magnetic field being distributed throughout the field. 
The analogy between Z, C and R in two-dimensional problems was introduced as a way of deducing 
these parameters from each other. 

The special case of energy storage in iron was discussed, leading to an expression for the work done 
in taking a sample of iron once around its hysteresis loop and to the idea of hysteresis loss. Finally, 
the consideration of the effects of a changing magnetic field on a solid conductor led to a discussion 
of eddy currents and of the losses associated with them. 
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7. Transmission lines 

Objectives 

• To introduce electromagnetic waves through a discussion of transmission lines as distributed 
circuits. 

• To discuss the solutions to the wave equation for sinusoidal waves and pulses, and to 
introduce the idea of phase velocity. 

• To explain the concept of characteristic impedance and its relationship to the power flow in 
the line. 

• To explain the use of complex notation to describe waves. 

• To show that a wave is reflected at the end of a line unless it is terminated by a load equal to 
its characteristic impedance. 

• To introduce the terms 'voltage reflection coefficient' and 'voltage standing wave ratio'. 

• To show how impedances are transformed by a transmission line, and to introduce the idea of 
the quarter-wave transformer as a matching device. 

• To consider the field description of a transmission line and to show how the phase velocity 
and characteristic impedance can be calculated. 

• To derive expressions for the electric and magnetic fields in a coaxial cable and to consider 
how the power flow in the line may be calculated from them. 

7.1 Introduction 

Two-wire transmission lines appeared early in the history of the practical application of electricity 
with the advent of the electric telegraph. They are still an important means of transmitting both 
electrical power and information, although for the latter purpose they are now supplemented by radio 
and microwave links and optical fibres. The commonest forms of transmission line are parallel- wire 
and twisted-pair lines, used for electricity distribution and telephone connections, and coaxial cables, 
used for television aerial downleads and for interconnecting electronic instruments. 

Throughout this book I have tried to show the links between the field and the circuit approaches to the 
description of electromagnetic phenomena. These provide alternative ways of dealing with problems. 
Electronic engineers are generally happiest when they are able to use circuit methods, with the 
properties of the components represented by equivalent circuits. Field theory provides an alternative 
to experimental measurements as a means of determining the equivalent circuit parameters. In this 
chapter, however, the field and circuit approaches draw even closer together. We shall see that the 
transmission of signals on transmission lines can be described in terms of propagation of 
electromagnetic waves. We shall also see that, when the dimensions of a circuit become comparable 
with the wavelength of these waves, elementary circuit theory breaks down. 



116 



Electromagnetism for Electronic Engineers 



7. Transmission lines 



7.2 The circuit tlieory of transmission lines 



It was shown in Chapters 2 and 6 that a two-wire transmission line such as a coaxial cable has shunt 
capacitance and series inductance uniformly distributed along its length. At low frequencies it is 
usually possible to ignore these impedances and treat circuits as sets of lumped components connected 
together by wires whose length has no effect upon the operation of the circuit. We must now 
investigate the circumstances in which this approximation is not valid and the consequences for the 
design of electronic circuits. 



V2Ldx V7L6X 

•^mnnnr — r^mRnnr o 



/ + d/ 



;cdx 



V ^ 6V 



Fig. 7.1 The 'tee' network representing a short length dx of a transmission line having series 
inductance L per unit length and shunt capacitance C per unit length. 



Figure 7.1 shows the circuit representation of an element of transmission line of length dx. If the line 
has inductance L per unit length and capacitance C per unit length, then the inductance and 
capacitance of the element must be L dx and C dx. 
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The inductance has been divided into two so that the network is a symmetrical tee section. For the 
moment we assume that the whole line stretches to infinity in either direction, being made up of 
identical tee sections joined together in a cascade. If the voltage across the capacitor is F, then 

r = V --L — dx = V + dV + -L — il + dl)dx 
2 dt 2 dr ^ 

didi) 

When the second-order term — ^ — - is neglected we obtain 

dt 

?^ = -L?L (7,) 

dx dt 

The current through the capacitor is 

dl = -Cdx^^ 
dt 

To first order V = V ,so that 

^ = -C^ (7.2) 

dx dt 

The pair of simultaneous differential equations (7.1) and (7.2) defines the relationship between the 
voltage and the current in the limit when the length (dx) of the element shrinks to zero and the line 
becomes completely uniform. 

The current can be eliminated from these equations by differentiating Equation (7.1) with respect to x, 
and Equation (7.2) with respect to t, since 

d^I d^I 



dxdt dtdx 



giving 



d V d^V 

^ = LC^- (7.3) 

dx' dt' 

Similarly, differentiating Equation (7.1) with respect to t and Equation (7.2) with respect to x, we get 
4 = ic4 (7.4) 
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Equations (7.3) and (7.4) are examples of the one-dimensional form of the wave equation. For 
sinusoidal waves their general solution has the form 

V = V^ cos{cot ± fcc) + ^2 sm{cot ± he) (7.5) 
and 

I = Ii co^{cot ± kx) + ^\n{cot ± kx) (7.6) 

where V\, F2, I\ and I2 are constants whose values are determined by the boundary conditions of the 
problem. Substituting these expressions into the wave equations shows that they are acceptable 
solutions provided that 

= co^LC {1.1) 

The physical significance of these solutions can be discussed, without loss of generality, by 
considering the behaviour of the function 

V = V^co^{cot-kx) (7.8) 

The form of this function is shown in Fig. 7.2. The solid curve shows the variation of voltage with 
position along the line when t = 0. At a later time t = t\ the wave has moved to the position shown by 
the broken line. This can be verified by noting that the first zero of the function is given by 

{(Dt - kx) = y 



X = 









yk J 




[ikj 



(7.9) 



This equation shows that the zero crossing, and hence the whole wave, is moving in the positive x- 
direction with a constant velocity given by 

This velocity is known as the phase velocity of the wave. 
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Fig. 7.2 The variation of the voltage along a transnnission line, as given by Equation (7.8) for f = 0 and 
a later tinne t = U showing how the wave nnoves steadily in the positive x-direction. 

From Equation (7.7) we also have 




(7.11) 



Equation (6.60) shows that Vp depends only on the properties of the material between the conductors 
and not on their geometry. 
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The propagation constant, or wavenumber, k is related to the wavelength X of the wave by 

271 

k = — (7.12) 

A 

as can be seen from Fig. 7.2. If k is negative the phase velocity is also negative and the wave is 
travelling in the negative x-direction. This corresponds to the positive signs in Equation (7.5). The 
alternative solution, involving the sine function in place of the cosine, gives a wave which is shifted in 
phase by 90^ 

By making a suitable choice of the constants V\ and V2 in Equation (7.5) we can make the amplitude 
and phase of the wave what we will. This is made clearer by a consideration of the alternative form of 
(7.5): 

V = VQCO^{cot±kx + (l)) (7.13) 

where the amplitude Fq and the phase (j) are arbitrary constants whose values are to be determined by 

the boundary conditions. The equivalence of the two forms of the wave given in Equations (7.5) and 
(7.13) can be demonstrated by expanding the cosine function in Equation (7.13). 

7.3 Representation of waves using connplex nunnbers 

In the theory of a.c. circuits 'j-notation' is an indispensable tool. Before proceeding to a discussion of 
the properties of transmission lines it is necessary to extend the notation to problems where the 
voltages and currents vary sinusoidally in distance as well as in time. To do this consider the voltage 

V = jV^sm{cot-kx)- jV^co^{cot-kx) (7.14) 

When this voltage is added to the general solution for the wave in the positive x-direction, obtained 
from Equation (7.5) by taking the negative signs, the result is 

V = v\co^{(ot-kx)+ j^\n{p)t-kx^- jV2\co^{cot-kx)+ j^m{(ot-kx^ 



or 

V = V exp[7(^^ - kx)] (7. 1 5) 

where the amplitude of the wave V is, in general, complex. It turns out that Equation (7.15) provides 
a more convenient form of representation for a wave than Equation (7.5), but it must be remembered 
that it is only the real part of Equation (7.15) which has physical significance. 
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Hence, whenever the form in Equation (7.15) is used to represent a wave in a problem it is implicit 
that, when an expression has been obtained as the solution to the problem, the real part of it is to be 
taken as having physical significance. It is easy to verify by substitution that Equation (7.15) is a 
solution of Equation (7.3). The two arbitrary constants in Equation (7.5) are present as the real and 
imaginary parts of V. 

7.4 Characteristic impedance 

So far we have discussed the propagation of waves on the line solely in terms of the voltage. The 
corresponding solution for the current can be obtained by substituting 



V = V^ Qxip[j{cot - kx)] 



and 



I = Ii Qxip[j{cot - kx)] 



into Equation (7.1), with the result 



-jkV,=-j(oLI^ 



so that, making use of (7.7) 



' k ^ 



(7.16) 



where 



2o = 



coL 



(7.17) 



This is a constant which has the dimensions of resistance. It is known as the characteristic 
impedance of the line. Equation (6.61) shows that Zq depends upon the geometry of the line and the 
properties of the material between the conductors. For the loss-less line considered here Zq is a real 
quantity, that is, a pure resistance, because the voltage and current are in phase with each other. It can 
be shown that if, on the other hand, the line has loss, then the voltage and current are no longer in 
phase with each other and, consequently, Zq has an inductive or capacitive component. The physical 
significance of the characteristic impedance becomes clearer when we consider the reflection of 
waves from the termination of the line. 
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7.5 Reflection of waves at the end of a line 



Figure 7.3 shows a transmission line of characteristic impedance Zq which is terminated by an 
impedance Z at x = 0. 

► Vf exp J(cot - kx) 



Vr exp j{(ot + he) 



Zo 



x = 0 

Fig. 7.3 A line of characteristic impedance Zq terminated by an impedance Z with incident and 
reflected waves having amplitudes \//and Vr. 

The wave incident on the termination is 



V = V.QiL^[j{cot-kx)\ 



(7.18) 



We must assume, until it has been shown otherwise, that some of the incident wave will be reflected 
as the wave 



(7.19) 
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The amplitudes Vi and Vr are complex. The corresponding currents are 



, 7 



QX^\j{(Ot - foe)] 



(7.20) 



and 



1 = 



^ 7 



exp[y ((^^/^ + Ax)] 



(7.21) 



as can be shown by substituting the expressions for the incident and reflected voltage waves in 
Equation (7.16). We shall see later that the negative sign in Equation (7.21) indicates that the power in 
the reflected wave is travelling in the negative x-direction. 

At X = 0 the voltage is 



F = (f^. + Fjexp(>0 



(7.22) 



and the current is 



/ = 



7 



exp(y^«;^) 



(7.23) 



by adding together the voltages and currents for the incident and reflected waves. 

But the ratio of the voltage to the current at x = 0 must equal the termination impedance Z. Therefore 



Z = ^ 



(7.24) 



which can be rearranged to give 



v.. Z-Z^ 



v.. z+z. 



<L=r 



(7.25) 



where F is known as the voltage reflection coefficient. Equation (7.25) shows that the amplitude of 
the reflected wave is zero when Z = Zq. When this condition is satisfied the termination is said to be 
matched to the line. 



When a line is terminated by a load equal to its characteristic impedance the instantaneous power 
absorbed in the termination is 
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W = VI = 



cos^ (cot) (7.26) 



The mean power is obtained by averaging Equation (7.26) over one complete cycle, giving 



i2 



W = ^ (7.27) 

Since there is no reflected wave when Z = Zq it follows that W is the mean power flow in the incident 
wave. This is the condition for maximum transfer of power from the incident wave to the load and it 
is a generalization of the Maximum Power Transfer Theorem of circuit theory. 

If the load resistor were replaced by a semi-infinite section of transmission line having the same 
characteristic impedance there would still be no reflected wave in the part of the line to the left of x = 
0. Thus the characteristic impedance is the input impedance of the semi-infinite line. Since we could 
just as well have replaced the load resistor by a section of line of finite length terminated by Zq 
without producing a reflected wave it follows that the input impedance of a transmission line 
terminated by its characteristic impedance is also Zq. 

When the line is not terminated by its characteristic impedance the voltage at a general point on it is 
V = V^ exp[7(^^ -kx)] + T V. exp[7 (^^ + kx)] (7.28) 
The amplitude of the voltage which would be measured at this point is 

I F| = I 1 1 exp(- jkx) 1 1 1 + r exp(2 jkx) \ = \V.\\l-\-r exp(2 jkx) \ (7.29) 
since | exp(- jkx) | = 1 . 

Equation (7.29) shows that the amplitude of the voltage varies along the line between the maximum 
and minimum values 

^■nax=|^|(l + |r|) (7.30) 
^.ax=|^|(l-|r|) (7.31) 

The ratio of these two is easily measured with a sliding probe and a suitable detector. It is known as 
the voltage standing wave ratio (VSWR), and its value is given by 
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1+ r 



s = 



(7 32) 



1- r 



The probe, therefore, detects a signal whose ampHtude varies periodically with position along the line 
as a standing wave. The wavelength of the standing wave observed is given by 2kx = 2n from 
Equation (7.29), that is 



where X is the wavelength of the travelling waves. 

7.6 Pulses on transmission lines 

So far we have assumed that the signals propagating on the transmission lines are sinusoidal. It is, of 
course, possible to consider non-sinusoidal signals by using Fourier synthesis to construct them from 
sine waves of different frequencies. But pulses, which are important for digital data transmission, can 
most easily be studied directly. It is easy to show by substitution that any functions having the form 



2 



(7.33) 




and 




(7.34) 



are solutions of Equations (7.1) to (7.4). 
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The function / can have any form, but for the study of 
pulses it is taken to be the step function shown in the figure 
on the right. This function is defined to be 0 for 
[x - Vpt) > 0 and 1 for (x - v^^) < 0 . Substituting the 

voltage and current given in Equation (7.34) into Equation 
(7. 1) gives 

VJ' = v^LIJ' 

where /' is the derivative of f so that, making use of Equations (7.11) and (7.17) 

^o=Vo (7-35) 
as before. 

Figure 7.4(a) shows a voltage source with internal impedance Zs connected to a semi-infinite 
transmission line of characteristic impedance Zq. If the line is uncharged until the switch is closed at 
t = 0 , then a wave having the form shown in Fig. 7.4(b) will propagate down the line with velocity Vp. 
As Vp is independent of frequency the step propagates without any change in its shape. A line for 
which this is true is said to be non-dispersive. Conversely, if Vp depends on frequency, the line is 
dispersive, the different Fourier components of the step travel at different speeds, and the shape of the 
step changes as it propagates. 

Zs 




(a) 

V 

VsZo/(Zo^Zs) 



(b) X 

Fig. 7.4 A semi-finite transmission line of characteristic impedance Zq is connected to a source of 
impedance Zs by a switch at f = 0. The result is a step voltage wave which travels down the 
transmission line as shown. 
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The voltage and current at the start of the hne must satisfy Equation (7.35), so that 

V,=-^V, (7 36) 

by the potential-divider rule. 



7.7 Reflection of pulses at the end of a line 

When a transmission line is terminated by a resistance R, part of an incident pulse is usually reflected. 
The voltage and current at any point on the line are given by 

(v-v) 

V = V, + V^ and 1 = 1+1 =^ d. 



I r 



To satisfy the boundary conditions at the termination 



V V.+V 

R = - = - — -z. 



SO that 



(7.37) 

From Equation (7.37) it is clear that there are three possible conditions: 

1 . R> Zo'.Vr is positive so that the reflected wave is added to the incident wave. 

2. R < Zq! is negative and the reflected wave is subtracted from the incident wave. 

3. R = Zq. Vr is zero and the incident wave is completely absorbed by the termination. 

Conditions (1) and (2) can be illustrated by the extreme cases of termination of a line by an open 
circuit and a short circuit. Figure 7.5 shows what happens when there is an open circuit. Figure 7.5(a) 
shows the situation just before the incident wave has reached the open circuit. When it reaches it a 
reflected wave is generated to satisfy the boundary conditions. From Equation (7.37) the amplitude of 
the reflected wave is equal to that of the incident wave. A short time later the situation is as shown in 
Fig. 7.5(b), with the reflected wave superimposed on the incident wave. 
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la) 



2V, 



(b> 



Open circuit 



Fig. 7.5 The reflection of a step wave by an open circuit. The wave is shown (a) just before it reaches 
the end of the line, and (b) just after it has been reflected. 
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In the same way Fig. 7.6(a) shows a wave approaching a short circuit. This time the wave is inverted 
on reflection so that the sum of the wave ampHtudes is always zero at the short circuit. A short time 
after the reflection of the wave the situation is as shown in Fig. 7.6(b). 



(a) 



(b) 



Short circuit 



Fig. 7.6 The reflection of a step wave by a short circuit. The wave is shown (a) just before it reaches 
the end of the line, and (b) just after it has been reflected. 



It is instructive to consider the effects of these reflections at the input of the transmission line. To 
make things simpler we will assume that the source is matched to the line so that Zs = Zq as shown in 
Fig. 7.7(a). Figure 7.7(b) shows how the voltage at the start of the line varies with time when the 
termination is an open circuit. Initially = /2 as the source feeds current into the line to set up the 
incident wave. This wave travels down the line in time T. At time 27 the reflected wave returns to A 
so that =V^ . The reflected wave is completely absorbed by the source impedance, so no further 
change takes place and the line is charged to the full source voltage. 




Fig. 7.7 Variation of the voltage at the input of a transmission line with time. The line has 
characteristic impedance Zq, and is connected to a matched source and to a mis-matched load. When 
the termination is an open circuit the voltage varies as shown in (b). (c) shows what happens when R 
is a short circuit. 
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If the line is terminated by a short circuit the voltage varies with time, as shown in Fig. 7.7(c). There 
is an initial period when the voltage at A is not zero, but this is cancelled at time IT hy the return of 
the reflected wave. According to elementary circuit theory the voltages would take their final values 
as soon as the switch was closed. The true situation can be described by saying that, until the return of 
the reflected wave, the source has no information about the magnitude of the terminating impedance. 
The step cannot travel faster than the speed of light, so it always takes a finite time for the system to 
reach a steady state. 

The transmission of a short pulse along a line can be investigated by considering the superposition of 
positive and negative step functions separated by a short interval of time. 

7.8 Transformation of innpedance along a transnnission line 

Consider a transmission line having characteristic impedance Zq which is terminated at x = 0 by an 
impedance Zl (see Fig. 7.8). At the point on the line which is / from the load, x = - 1 and the voltage 
is given by 




(7.38) 



The corresponding expression for the current is 




(7.39) 



where 



0 



0 



(7.40) 



from Equation (7.25). 



The apparent impedance at this point on the line is given by 



^ exp(y^/)+rexp(-y^/) ^ 
^ / exp(y^/)-rexp(-y^/) ^ 



0 



(7.41) 
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After substituting for F and rearranging, we get 




[exp(#/ ) + exp(- jkl )] + Zq [exp(yl/ ) - exp(- jkl )] 
[exp(7l/ ) - exp(- jkl )] + Z^ [exp(7l/ ) + exp(- jkl )] 
Z^ cos(^/)+ yZ^ sin(^/) 



(7.42) 



Zq cos(^/)+ yZ^ sin(^/) 
Z^+yZp tan(^/) 

Zq + yZ^ tan(^/) 



This is a very important result. It shows that if the length / of the transmission line is of the same order 
of magnitude as the wavelength {Inlk), or greater, then the apparent impedance at the input of the line 
is not normally equal to the load impedance. The usual assumption of circuit theory is that the 
impedances are independent of the lengths of the connecting wires. We now see that this is just the 
limit of Equation (7.42) when kl tends to zero. 

The significance of Equation (7.42) can be explained in a slightly different way by considering Fig. 
7.8. This shows two transmission lines, one of which is longer than the other by /. The longer line is 
terminated by a load Zl and the shorter by Z^ , where Z[ is given by Equation (7.42). At the plane A-A 
the impedance presented to the incident wave is the same on both lines. Thus the reflected wave to the 
left of A-A must have the same amplitude and phase on both lines. It follows that Zl andZ^ must 

reflect the incident wave with the same amplitude but with a phase difference of exp (2jkl), since that 
is the phase change from A-A to Zl and back. 
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Fig. 7.8 The apparent ternninating innpedance of a transnnission line varies with position along the line. 
If two lines are ternninated by Z/. and Z'i_ as calculated fronn Equation (7.42), then the annplitude and 
phase of the reflected wave to the left of the plane A-A is the sanne on both lines. 

7.9 The quarter-wave transformer 

A case of considerable practical importance is obtained when the line is exactly one quarter of a 
wavelength long. Setting^/ = ;r/2 in Equation (7.42) gives 



7 7 



or 



Zl=Z,Z', (7.43) 

This result allows us to match a load to a line by inserting a short section of a line, having a different 
impedance, between them, as shown in Fig. 7.1 1. The short length of line must be a quarter 
wavelength long, and it must have the characteristic impedance given by (7.43). The same technique 
can be used to ensure a match between two transmission lines which have different impedances. 



2r 




V A 

Fig. 7.9 A source of impedance Z'/. can be matched to a load of impedance Z/. by using a quarter 
wavelength of a transmission line whose characteristic impedance is Zq = V(Z/. Z'/.) as a transformer. 
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The quarter- wave transformer suffers from the difficuhy that it works exactly only at a single 
frequency. This problem can be overcome by cascading a number of quarter-wave transformers to 
produce a broad-band match. 

Two special cases arise when the terminating impedance of the transformer is either an open circuit or 
a short circuit. If there is an open circuit, then = go and (7.42) becomes 



Zq ytan(^/) 

so that Ti tends to zero as kl tends to nil. Thus a quarter-wave transformer transforms an open circuit 
into a short circuit. The converse is also true. If the termination of the transformer is a short circuit, 
then the impedance at its input tends to infinity. In either case there is a total reflection of the incident 
power. 



The simplest practical transmission line from the field point of view is the coaxial cable. We shall 
therefore take it as the example to be discussed, but the methods applied, and the conclusions reached, 
are applicable to all types of uniform two-wire line. A typical cable has a solid copper inner conductor 
and a braided outer conductor separated by a uniform cylindrical layer of polythene. To make the 
problem tractable we assume that the braided outer conductor can be represented with sufficient 
accuracy by a continuous conducting cylinder. Taking the radii of the inner and outer conductors as a 
and h, respectively, we recall that the capacitance per unit length is given by (6.50) 




1 



7.10 The coaxial line 



Ins 



(7.44) 




and the inductance per unit length is given by (6.37) 




(7.45) 



Substituting these expressions for C and L into Equations (7.11) and (7.17) we get 



1 



(7.46) 




and 
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(7.47) 



A very interesting result is obtained for the special case of an air-spaced line when 



= (8.854xl0"''x4;rxl0"')"2 
= 0.2998x10' ms"' 



which is the experimental value for the velocity of light in free space. This suggests the existence of a 
link between electromagnetism and optics, and the possibility that light is an electromagnetic 
phenomenon. Since the relative permittivity of any medium is greater than unity, it follows that the 
phase velocity of the waves on a transmission line is always less than or equal to the velocity of light. 

7.11 The electric and magnetic fields in a coaxial line 

To find expressions for the electric and magnetic fields within a coaxial line which is carrying a wave 
along it we apply Gauss' theorem and the magnetic circuit law. This approach involves us in making some 
assumptions about the fields in the line. Taking the electric field first, we assume that, as in the electrostatic 
case, only the radial component is present. The charge per unit length on the centre conductor is given by 




(7.48) 
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so, applying Gauss' theorem 



q _CV,\ 



Insr 2ns r 



exp[y((^^/^-Az)] 



Substituting for the capacitance per unit length from Equation (7.44) gives 



= , ^ — exp[y ((^^ t-kz)] 



hi{b/a)r 



(7.49) 



The magnetic field is given by 



1 F„ 



exp[7(<y?-A:z)] 



Substituting for Zq from Equation (7.47) 



(7.50) 



These fields are shown in Fig. 7.10. We note that the ratio of the strengths of the electric and magnetic 
fields is 



= Z,. 



(7.51) 



This quantity, which depends only upon the properties of the material filling the line, has the 
dimensions of resistance. It is known as the wave impedance of the wave. 




Electric fietd [E) 
Magnetic fieJd jw^^ 
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Fig. 7.10 The electronnagnetic field in one wavelength of a coaxial line showing the directions of the 
electric and nnagnetic fields and the direction of propagation of the wave. Half of the outer conductor 
has been rennoved for clarity. The electric field lines are shown schennatically; they actually radiate 
equally in all directions. 

Although these results have been derived by making assumptions about the directions of the fields, 
they are confirmed by more rigorous analysis. It is found that, for any two-wire transmission line, the 
electric and magnetic field vectors are perpendicular to each other and to the direction of propagation 
of the wave. Such waves are known as transverse electric and magnetic (TEM) waves. 

7.12 Power flow in a coaxial line 

When discussing the storage of energy in capacitors and inductors we saw that it is sometimes useful 
to think of the energy as being distributed throughout the electric or magnetic fields. In the same way 
we can think of the power flowing down a coaxial cable as being distributed throughout the field. At 
any point the electric energy density is 

1 1 9 

=-B'E = -sE^ (7.52) 
while the magnetic energy density is 

^M=\s-H = ^jU,Hl (7.53) 
Making use of Equation (7.51) we can write the total energy density as 

= £E^ (7.54) 



1 1 ^ - ^ 



2 ' 2 



\Mo J 



Now at any point within the line the fields vary sinusoidally with time. Therefore the time average of 
the stored energy is 

w = -sE]= , (7-55) 
2 ' 2(ln(Va)) 

The whole field pattern is moving in the positive z direction with velocity v^, so the power density is 
S = wv^ J-eXsv,e)=]^E^H, (7.56) 
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making use of Equations (7.46) and (7.51). We shall see in Chapter 8 that this is an example of the 
application of a more general result relating to the power flow in electromagnetic waves. 

To check the equivalence of the field and circuit descriptions of the problem we can calculate the total 
power flow by integrating Equation (7.56) over the space between the conductors. Then 



from the expression for Zq given in Equation (7.47), so the power is identical to that calculated by 
circuit methods (Equation (7.27)). 

7.13 Summary 

In this chapter transmission lines have been introduced through a discussion of distributed circuits. 
The voltages and currents were found to be governed by the wave equation, and it was demonstrated 
that solutions could be found which could be interpreted as travelling waves. The power carried by 
these waves was expressed in terms of the voltage on the line and its characteristic impedance. 




dr 



(7.57) 




(7.58) 
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It was demonstrated that, unless a line is matched by being terminated by a load equal to its 
characteristic impedance, some of the power in a wave incident upon the termination is reflected back 
down the line. This reflection produces a partial standing wave on the line. The effect of the reflected 
wave is to make the input impedance of a line depend upon its length as well as upon the terminating 
impedance. Thus it was shown that, at high frequencies when the lengths of lines are comparable with 
the wavelengths of the signals on them, it is no longer possible to treat circuits in terms of lumped 
components connected together by wires of arbitrary length. The transformation of impedance by a 
short length of a transmission line is put to use in the quarter wave transformer to match a load to a 
line. 

It is also important for lines to be correctly matched when the signals are in the form of short pulses. 
Incorrect matching results in reflections of the pulses. If a line is not matched at both ends multiple 
reflections can occur so that a whole train of pulses is produced by transmitting a single pulse down 
the line. 

Transmission lines were also considered from the point of view of the electric and magnetic fields 
within them. It was shown that the phase velocity and characteristic impedance of a line at a given 
frequency can be calculated from field considerations. Expressions were found for the electric and 
magnetic fields within a coaxial line, and it was demonstrated that the flow of power down the line 
could be regarded as that of an electromagnetic wave propagating in the space between the 
conductors. It was found that, for an air-spaced line, this wave propagates with a phase velocity equal 
to the speed of light. The equivalence between the circuit and field approaches to the description of 
waves on transmission lines was demonstrated by showing that both give the same expression for the 
power flow. 
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8. Maxwell's equations and electromagnetic 
waves 

Objectives 

• To show how Maxwell removed an inconsistency in the magnetic circuit law by introducing 
the idea of the displacement current. 

• To derive the differential forms of the magnetic circuit law and Faraday's law, and to 
introduce the curl of a vector. 

• To present Maxwell's equations in both integral and differential form. 

• To demonstrate how the existence of plane electromagnetic waves in free space can be 
deduced from Maxwell's equations. 

• To show that the power flow in an electromagnetic field can be calculated by integrating the 
Poynting vector over a closed surface. 

8.1 Introduction 

In Chapter 7 we saw how it is possible for electromagnetic waves to propagate along a transmission 
line such as a coaxial cable. We shall now proceed to show that electromagnetic waves can exist 
independently of any system of conductors and can travel through free space. This conclusion was 
first published by James Clerk Maxwell in 1873. The set of equations which bears his name forms a 
summary of all the topics treated in this book and demonstrates the symmetry and unity of 
electromagnetic theory. At the same time they are the foundation upon which modem electromagnetic 
theory is built. We are now so accustomed to the presence of radio, television, mobile phones and 
radar in the world that it is hard to realize how great Maxwell's contribution was to the understanding 
of the theory of electromagnetism. In his own day his theory was regarded as quite outrageous by 
many until its validity was demonstrated by Hertz's experiments some fifteen years later. 

8.2 Maxwell's fornn of the nnagnetic circuit la 

Maxwell noticed that, under certain circumstances, the magnetic c 
law in the form derived in Chapter 4 can give inconsistent results, 
is readily demonstrated by considering the circuit shown in the fig 
the right - a charged capacitor discharging though a resistor. 



To apply the magnetic circuit law we must define a closed path C 
which encircles the wire and an open surface S which spans that 
path as a soap film can span a loop of wire. The figure in the 
margin attempts to show what is meant. If S is chosen so that it 
passes through the wire, then no problem arises and the magnetic 
field satisfies Equation (5.7), which is repeated here for 
convenience: 
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|Hdl = jJjdA 



(8.1) 



c s 



But, if the surface S is chosen so that it passes between the plates of the capacitor, then J = 0, 
implying that the line integral of H around C is also zero. This sounds nonsense! We expect the 
magnetic field around the circuit to be independent of the details of the means by which it has been 
calculated. 

To resolve this difficulty we consider two different surfaces Simd S2 which both span C, as shown in 
Fig. 8.1, and which together form a closed surface. 



Fig. 8.1 Maxwell showed that the conduction current (/) passing through a closed loop (C) can be 
dependent upon the choice of the surface (S) over which it is calculated. The sum of the conduction 
and displacement currents (/ + d^/dt) is independent of the choice of surface. 
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At a particular instant let the total charge enclosed within this surface be Q, and let currents Ii and I2 
cross the two surfaces as shown. Now charge must be conserved, so the rate of change with time of 
the charge enclosed must be equal to the net current flow into the volume enclosed. That is, 

/,-/,=f (8.2) 

at 

If we define Wi to be the flux of D through 5*1 and W2 that through 5'2,as shown in Fig. 8.1, then the 
application of Gauss' theorem gives 

^i-^2=e (8.3) 

Differentiating Equation (8.3) with respect to time and substituting for dQ/dt in Equation (8.2) gives 
^ d^, ^ d^. 

This equation shows that the quantity (/ + d^/dt) is independent of the choice of the position of the 
surface S even when / is not. This quantity is known as the total current, and d^/ dt is known as the 
displacement current. 

When the current density on the right-hand side of Equation (8.1) is replaced by the total current 
density (j + dD/dt), the result is the equation 

lH-dl=fffj + — 1-dA (8.5) 
J ^ dt J 

This form of the magnetic circuit law is universally true. 

8.3 The differential form of tine magnetic circuit law 

In Chapter 1 we saw that it is possible to express Gauss' theorem in either an integral form (Equation 
(1.5)) or a differential form (Equation (1.9)). The differential form turns out to be the more useful of 
the two because it enables a wider range of problems to be solved. The same is true of the magnetic 
circuit law, as we shall see. 

To derive the differential form of Equation (8.5) we apply it to the small rectangular path ABCD 
shown in Fig. 8.2. The surface S is taken to be in the j^-z plane. The line integral of H is taken around 
ABCD in a direction which is in the right-handed corkscrew sense with respect to the positive x- 
direction. In other words, the elementary loop is encircled in the direction ^5 CD. 
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Fig. 8.2 The differential fornn of the nnagnetic circuit law can be derived by considering the line integral 
of the nnagnetic field around the infinitesinnal loop ABCD. 



The value of the line integral of H around the loop is computed by considering each side in turn. The 
effect of the dot product of H with dl is to pick out the component of H which is parallel to the path. 
Thus, for AB we get 



(8.6) 



Along CD the j;-component of H can be written {h^ + {dH^ j dz)&) , so that 



J^H.dl = 



dH 

f 

dz 



& 



(8.7) 



The minus sign on the right-hand side of this equation is there because the path CD is traversed in the 
negative j^-direction. 

Combining Equations (8.6) and (8.7) gives 




(8.8) 



Applying the same argument to the other two sides of the loop gives 



\ H dl + f H dl = — '-dydz 

Jb Jd Ph, 



dy 



(8.9) 



The complete line integral around the loop ABCD is found by adding Equations (8.8) and (8.9) to give 



H dl 



dy dz 



dy dz 



(8.10) 
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Since the loop is small, it is reasonable to assume that J and D are constant over it, so that the integral 
on the right-hand side of Equation (8.5) becomes 



ni 



, dD 

dt 



dA = 



dt 



dy dz 



(8.11) 



This time the dot product picks out the components of J and D which are normal to the plane of the 
loop because the vector area dA is also normal to the plane of the loop. 

Equating Equations (8.10) and (8.1 1) yields the differential form of Equation (8.5) for the x-direction: 



■ = J+- ^ 



dy dz 



dt 



(8.12) 



Similar equations can be obtained for the other two coordinate directions: 



dH^ dH^ , dD^ 

i =- = J+ — ^ 

dz dx dt 



(8.13) 



dH^ dH^ ^ dD^ 



dx dy 



dt 



(8.14) 
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Equations (8.12) to (8.14) together form the differential version of Equation (8.5) expressed in terms 
of the components of the vectors. 

As usual, it is possible to write these equations in a more compact form by using vector notation. To 
do this we multiply each by the appropriate unit vector and add them together to give 



dy dz 



x + 



dz dx 



y + 



dH^ dH^ 



dx dy 



z = J + — 

dt 



(8.15) 



The left-hand side of this equation is known as the curl of H and can be written as a determinant: 



X 


y 


z 


d 


d 


d 


dx 


dy 


dz 









curlH=— — — (8.16) 



which can be recognized as the vector product of the vector operator V and the vector H. Thus 
Equation (8.15) can be written succinctly in the form 

VaH = J + — (8.17) 
dt 

Like the other differential forms derived earlier, Equation (8.17) is actually valid for coordinate 
systems other than the rectangular Cartesian system used to derive It. 

8.4 The differential form of Faraday's law 

The general integral form of Faraday's law was shown to be 



E.dl = -ff^dA 

JJ dt 



(8.18) 



(Equation (6.18)) for the case when the circuit is fixed in space and the magnetic flux density is 
changing with time. 



So far we have assumed that the line integral of E is to be taken around a loop of wire. But if the loop is 
open-circuited and the potential difference between its ends is measured with a. high-impedance volt- 
meter, then the loop can be thought of as a device for measuring the value of the line integral of an electric 
field which exists in space whether the wire is present or not. Thus Equation (8.18) can be regarded as a 
generalization of Faraday's law, implying that, if there is a changing magnetic field in any region of space, 
then there is also an electric field there. This is, strictly speaking, a plausible guess rather than a deduction, 
but it is one whose validity has been demonstrated by the correctness of the results derived from it. 
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Making a comparison between Equations (8.5) and (8.18) allows us to deduce straight away that the 
differential form of (8.18) must be 

VaE = (8.19) 

dt 

8.5 Maxwell's equations 

We have now introduced all the principles of electromagnetism and it is convenient to gather the main 
results together. Taking the integral forms first they are: 

Gauss' theorem in electrostatics (Equation (2.5)) 

(^D-dS = (8.20) 
Gauss' theorem in magnetostatics (Equation (4.15)) 

(^B-dA = 0 (8.21) 
The magnetic circuit law (Equation (8.5)) 

1 H.dl = jj|^J + ^j.dA (8.22) 
Faraday's law of induction (Equation (6.18)) 

fE-dl = -n^-dA (8.23) 
The corresponding differential forms are 

(2.6) V-D = p (8.24) 

(4.16) V-B = 0 (8.25) 

dD 

(8.17) VaH = J + (8.26) 

dt 

8R 

(8.19) VaE = (8.27) 

dt 
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These four equations form a summary of the whole of fundamental electro-magnetic theory. They are 
known, collectively, as Maxwell's equations and are the starting point for the discussion of all the 
more advanced topics in electro-magnetism. 

In order to make use of them we also need a number of other equations which have been introduced in 
previous chapters. First there is the continuity equation 



(3.14) V-J = -i^ (8.28) 

dt 

Then there are the constitutive relations which introduce the properties of materials, albeit in an 
idealized form. These are: 

(2.4) D = ^E (8.29) 
(3.3) J = crE (8.30) 

(5.5) B = //H (8.31) 



The application of Equations (8.24) to (8.31) is most easily accomplished by making use of the 
methods of vector calculus. These techniques lie beyond the scope of this book, but a simple example 
of the use of Maxwell's equations which can be dealt with using more elementary mathematical 
methods is the subject of the next section. 
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8.6 Plane electromagnetic waves in free space 

In the previous chapter we saw that electromagnetic waves can propagate along a coaxial line and that 
the electric field vector, the magnetic field vector, and the direction of propagation are all mutually 
perpendicular. Taking our clue fi-om this, we can see whether Maxwell's equations indicate that 
similar waves can propagate in free space far from any material boundaries. 

Let us assume that the wave propagates in the z-direction, and that the electric field has a component 
only in the x-direction. We assume, furthermore, that the intensity of the electric field varies only in 
the z-direction. Then from Equation (8.27), making use of the definition of the curl of a vector given 
in Equation (8.16), we have 

y^ = -^ (8.32) 

^ dz dt 



since all the other components of V x E are zero. This equation shows that the magnetic field vector 
must lie in the j-direction, at right angles to both the electric field and the direction of propagation. 
Thus Equation (8.32) can be written as a scalar equation 

— = (8.33) 

dz dt 

In free space the conduction current J = 0, so Equation (8.26) becomes 
dH^ dD 

-^ = ^ (8.34) 

dz dt 

The magnetic field terms can be eliminated between these two equations by differentiating the first 
with respect to z and the second with respect to t, and by making use of Equations (8.29) and (8.31) 
with s = Sq and jU = /Uq. The result is that 

That is 

d'E 1 d'E 

^ = ^ f- (8.36) 

dz' v\ dt' 



where 
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v^=^= (8.37) 



is the velocity of light. Equation (8.36) is the wave equation which was first encountered in Equation 
(7.4). 

If the electric field terms are eliminated from Equations (8.33) and (8.34) in a similar manner, the 
result is 

d'H^ 1 d'H^ 

f = ^ (8-38) 

The two wave equations have the following general solutions for waves travelling in the positive z- 
direction: 

E^=E,^^^[j{(Dt-kz)] (8.39) 
and 

H^=H,Qxv[j{cot-kz)\ (8.40) 

where and Hq are complex amplitudes which incorporate the relative phases of the electric and 
magnetic fields. 

The relationship between the two fields can be found by substituting Equations (8.39) and (8.40) into 
either Equation (8.33) or Equation (8.44). Making use of Equation (8.33), we get 

- jkE^ exp [j{cot - kz)] = -jcOjU.H, exp [j{cot - kz)] (8.41) 
From this it can be seen that and Hy are in phase with one another. Also that 

^0 _ I Mo 



Ho 



= 377 n (8.42) 



We recognise this as the wave impedance of a wave in free space (see Equation (7.51)). It is 
sometimes referred to as the intrinsic impedance of free space. 

Figure 8.3 illustrates the complete solution which has been obtained. It is necessary to be a little 
careful in interpreting this diagram. It shows the amplitudes of the field vectors rather than being a 
map of the field. The whole field pattern is moving in the z-direction with the velocity of light. At any 
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instant the electric and magnetic field strengths are uniform at all points in a plane perpendicular to 
the z-axis. 




Fig. 8.3 The electric and magnetic fields in a plane electromagnetic wave propagating in the z- 
direction. 



150 



Electromagnetism for Electronic Engineers 



8. Maxwell's equations and electronnagnetic waves 



8.7 Power flow in an electromagnetic wave 

In Chapter 7 it was shown that the power flow in a coaxial Hne could be calculated correctly by 
assuming that the stored energy in the electromagnetic field travelled with the phase velocity in the 
direction of the wave. For waves in free space the energy density is, from Equation (7.54), 

w = s,El (8.43) 

The instantaneous power density in the wave is 

S = wv^=^^ = E^H^ (8.44) 



Now S is in the z-direction, so it is possible to write Equation (8.44) in vector notation as 

S=EaH (8.45) 

The vector S is known as the Poynting vector. Poynting's theorem states: 

The integral of (E a H) over a closed surface is equal to the instantaneous flow of electromagnetic 
power out of the volume enclosed by that surface. 

The proof of this theorem is beyond the scope of this book. Note that, strictly speaking, it is the 
integral of S over a closed surface which has meaning rather than S itself We have already seen that 
the integral of the Poynting vector gives the correct answer for the power flow in a lossless coaxial 
line. 

8.9 Sunnnnary 

We have considered how Maxwell removed an inconsistency in the magnetic circuit law by 
introducing the idea of the displacement current. The modified law was expressed in differential form 
by introducing the curl of the vector H. Faraday's law was also put into differential form by using the 
curl of E. 

The mathematical statements of the laws of electromagnetism were collected together in both their 
integral and differential forms. This set of equations is known as Maxwell's equations. It was 
demonstrated that, for the special case of plane waves in free space. Maxwell's equations lead 
naturally to the prediction of the existence of electromagnetic waves which travel with the speed of 
light. 

Finally, it was shown that the power flow in the electromagnetic field can be represented by the 
Poynting vector (E a H) . 
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Physical constants 



Primary electric constant (80) 


8.854 X 10-^2 


Fm"' 


Primary magnetic constant (no) 


431 X 10"^ 




Hm"' 


Velocity of light in vacuum 


0.2998 X 


10^ 


m s' 


Wave impedance of free space 


376.7 




Q 


Charge on the electron 


-1.602 X 


10-19 


C 


Rest mass of the electron 


9.108 X 10"^' 


kg 


Charge/mass ratio of the electron 


-1.759 X 


10" 


Ckg- 



Properties of dielectric materials 

Relative permittivity 



Alumina 99.5% 


10 


Alumina 96% 


9 


Barium titanate 


1200 


Beryllia 


6.6 


Epoxy resin 


3.5 


Ferrites 


13-16 


Fused quartz 


3.8 


GaAs (high resistivity) 


13 


Nylon 


3.1 


JTCllCllllll W Cl A 




Perspex 


2.6 


Polystrene 


2.54 


Polystyrene foam 


1.05 


Polythene 


2.25 


PTFE (Teflon) 


2.08 


Properties of conductors 






Conductivity (S m"') 


Aluminium 


3.5 X 10^ 


Brass 


1.1 X 10^ 


Copper 


5.7 X 10^ 


Distilled water 


2 X 10"^ 


Ferrite (typical) 


10-^ 


Fresh water 


10-^ 


Gold 


4.1 X 10^ 


Iron 


0.97 X 10^ 


Nickel 


1.28 X 10^ 


Sea water 


4 


Silver 


6.1 X 10^ 
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Steel 0.57 x 10^ 

Properties of ferromagnetic materials 

|ir Saturation magnetism (Bsat) (T) 

FeroxcubeS 1500 0.2 

Mild steel 2000 1.4 

Mumetal 80 000 0.8 

Nickel 600 

Silicon iron 7 000 1.3 

Summary of vector formulae in Cartesian coordinates 

a A b = (a^b^ - a^b^ )x + {a^b^ - ajj^ )y + [ajj^ - a^b^ )z 

^ d ^ d ^ d ^ 
V = — x + — y + — z 

dx dy dz 

grad V = VV = x + y + z 

dx dy dz 
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^ da^ da da^ 
div a = V • a = — - + — - + - 



dx dy dz 

curl a = V A a = <^ — \x + \ — \y + <^ — \z 

ydy dz j [ dz dx } [ dx dy j 

. eV aV 

V V = — --\- — --\- — ^ 
dx' dy' dz' 



Summary of the principal formulae of electromagnetism 

Inverse square law of electrostatic force 

Relationship between E and D 

D = = s^s;e 

Gauss' theorem 

D.dS = jjj pdv divD = V.D = /7 

Electrostatic potential difference 

V.-V.=-{^ E'dl E = -gradF = -VF 

Poisson's equation 



72j, _ P 



^0 

Energy density in an electric field 

w = -D-E 

2 

Ohm's law 

J = crE E = pJ 

Power dissipated per unit volume 
P = E J 

Continuity equation 

Electromotive force 
^=1 E dl 

Law of force between moving charges 



F 



= ^%r + ^^^(v,A(v,Ar)) 



Force on a moving charge 
F = e(E + (vAB)) 
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Force on a current-carrying conductor 
F = /(dlAB) 
Biot-Savart law 
g /V[r dlAT 
471 ^ 

Relationship between B and H 
B = //H = //o//,H 

Magnetic circuit law as modified by Maxwell 



H.dl = JJ J + 



dt 



dA curlH = VAH = J + — 

dt 



C 5 V 

Conservation of magnetic flux 

<^BdA = 0 divB = VB = 0 

Faraday's law of induction 

f f r OB 5B 

I E-dI = -rr — -dA curlE = VAE = -— 

J J-' dt dt 

Energy density in a magnetic field 

w = -B H 

2 

The Poynting vector 
S = EaH 
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